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Preface 


This  report  summarizes  the  work  carried  out  from  7/15/82  -  2/28/90  under  ONR 
sponsorshipbn  laser-stimulable  luminescence.  The  effort  emphasized  three  aspects: 

1 .  Investigation  of  basic  laser-stimulated  luminescence  phenomena  such 
as  laser-heated  thermoluminescence  for  radiation  dosimetry  and  optical 
(non-thermal)  stimulated  luminescence  for  potential  application  in 
radiation  imaging  and  personnel  dosimetry. 

2.  Graduate  student  training  and  interaction  with  guest  scientists. 

3.  Technology  transfer  to  U.S.  industry. 

Dr.  Herschel  S.  Pilloff  was  the  technical  monitor.  We  are  greatly  indebted  to  him  for  his 
support  and  enthusiastic  interest  in  this  work,  and  particularly  for  his  encouragement  to 
apply  the  gained  knowledge  toward  solutions  of  problems  in  radiation  protection  in  U.S. 
industry  and  the  DOD. 


Introduction 


This  research  program  on  laser-stimulated  luminescence  was  conceived,  planned, 
and  executed  with  the  intent  to  explore  the  basic  physical  principles  of  laser-stimulation  of 
storage  phosphors  (those  that  retain  a  small  fraction  of  the  energy  absorbed  when  exposed 
to  ionizing  radiation)  and  their  applications  in  radiation  dosimetry  and  radiography.  From 
the  outset,  cross-fertilization  between  this  research  project  and  a  simultaneous  engineering 
development  program,  carried  out  at  International  Sensor  Technology,  Inc.,  (1ST,  Inc.) 
was  an  important  aspect  of  the  effort.  Information  on  the  physics  of  laser  heating  of  two- 
layer  structures  (the  phosphor  layer  on  a  substrate)  influenced  the  design  of  equipment  and 
thermoluminescence  dosimeter  configurations  at  1ST,  Inc.  which,  in  turn,  supplied  to  the 
ONR-sponsored  program  at  Washington  State  University  electronic  circuit  and  electro- 
optical  designs  for  its  experimental  setups.  This  collaboration  between  an  academic 
institution  and  U.S.  industry  has  been  sanctioned,  encouraged,  and  supported  by  ONR  as 
it  became  clear  that  the  work  would  yield  the  scientific  basis  for  the  development  of  a  new 
generation  of  instruments  and  dosimeters  for  personnel  dosimetry  to  be  used  by  the  U.S. 
Navy.  The  success  of  this  approach  has  been  extraordinary:  1ST,  Inc.  has  just  concluded 
an  Advanced  Development  Contract  for  the  U.S.  Navy  (Contract  No.  N60921-88-C-0085, 
NSWC,  Silver  Spring,  MD),  and  a  full  scale  engineering  contract  will  be  awarded  to  the 
successful  bidder  in  September  of  1990.  Thus,  all  prerequisites  for  use  of  the  new 
dosimetry  equipment  by  the  U.S.  Navy  are  in  place,  and  a  purely  academic  research  project 
will  have  been  transferred  to  industry  for  production  for  instrumentation  that  is  of 
importance  for  our  national  defense. 

Finally  we  mention  here  two  additional  spin-offs:  a  two-dimensional  dose  mapping 
apparatus  and  a  fiber-optic  dosimetry  probe.  Both  benefited  greatly  from  the  ONR 
supported  research  on  laser-heated  thermoluminescence.  These  instruments  are  presently 
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under  development  at  1ST,  Inc.  with  Small  Business  Innovation  and  Research  (SBIR) 
support  from  the  National  Cancer  Institute. 

A  chronology  of  events  in  the  development  of  all  technology  that  is  based  on  laser¬ 
heating  of  thermoluminescent  phosphors  is  given  in  Table  1. 

After  completion  of  the  work  on  laser-heated  thermoluminescence  we  investigated 
non-thermal  laser  stimulation  of  inorganic  phosphors  that  we  considered  to  be  potential 
candidates  for  dosimetric  applicadons.  This  part  of  the  effort  concentrated  on  two 
materials:  ZnS:Dy  for  use  in  neutron  radiography  and  CaS04:Sm  as  the  first  non-thermallv 
laser-srimulable  phosphor  for  personnel  dosimetry. 

The  work  on  ZnS:Dy  was  disconrinued  after  it  became  clear  that  neither  WSU  nor 
1ST,  Inc.  were  equipped  for  experiments  on  neutron  radiography  and  no  other  U.S.  firm 
indicated  any  interest  in  the  developed  concept.  A  patent  disclosure,  submitted  to  the  WSU 
Patent  and  Trademark  Office  was,  therefore,  not  acted  upon. 

The  development  of  the  first  laser-stimulable  phosphor  with  genuine  potential  for 
applications  in  personnel  dosimetry  must  be  considered  the  basis  for  future  work.  As  its 
laser  stimulation  is  non-thermal,  several  disadvantages  of  laser-heated  thermoluminescence 
may  be  eliminated:  heating  is  relatively  slow  as  compared  to  optical  stimulation  and  both, 
binder  and  substrate,  of  a  thermoluminescent  dosimeter  element  must  withstand  300- 
400°C.  An  optically  stimulable  dosimetric  phosphor  allows  the  construction  of  dosimeters 
using  organic  binders  for  the  sensitive  phosphor  layer  as  well  as  polymer  substrates, 
resulting  in  greater  simplicity  and  reduced  production  costs. 

In  the  following  we  first  present  the  main  results  of  the  laser-heating  effort  in 
thermoluminescence  dosimetry  (in  form  of  reprints  of  four  consecutive  publications)  and  a 
brief  summary  of  the  CaSO^Sm  work  together  with  a  copy  of  a  Master's  degree  thesis  on 
this  topic. 
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Part  I:  Laser-Heated  Thermoluminescence 


This  part  of  the  research  program  commenced  with  a  detailed  study  of  heating  a 
single  LiF  dosimeter  chip  and  a  two-layer  dosimeter  configuration  with  a  TEMoo  mode 
(Gaussian  beam  profile)  CO2  laser  beam  and  the  corresponding  thermoluminescence 
emission.  The  main  results  are  reviewed  in  the  paper  entitled  "Laser-Stimulated 
Thermoluminescence,"  presented  below. 

For  true  dosimetric  applications,  the  laser  beam  must  be  of  uniform  intensity 
profile.  Using  a  technique  invented  at  1ST,  Inc.,  we  continued  the  investigation  of  TL 
heating  with  a  uniform  square  beam.  This  work  is  reviewed  in  the  publications  following 
the  one  on  Gaussian  beam  heating. 
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Laser  stimulated  thermoluminescence 

A.  Abtahi  and  P.  Br&unlich 

Department  of  Physics,  Washington  State  University,  Pullman,  Washington  99164-2814 

P.  Kelly 

Division  of  Microstructural  Sciences,  Motional  Research  Council.  Ottawa,  Canada  K1A  0R6 
J.  Gasiot 

C.EA.'.,  Universite  des  Science  et  Techniques  du  Languedoc.  34060  Montpellier  Cedex,  France 
(Received  2  January  1985;  accepted  for  publication  28  March  1985) 

Experimental  and  computational  methods  are  presented  for  the  complete  characterization  of  the 
thermoluminescence  response  obtained  from  thermoluminescent  phosphors  upon  exposure  to 
localized  Gaussian  laser  heating  beams.  A  number  of  different  phosphor  configurations  are 
described  as  examples.  These  include  LiF:Mg,Ti  (TLD-100,  Harshaw  Chemical  Corporation)  in 
form  of  chips,  which  are  widely  used  in  the  dosimetry  of  ionizing  radiation,  and  thin-layer 
dosimeters  prepared  either  as  self-supporting  films  or  powder  in  a  polyimide  matrix,  or  on 
substrates  of  LiF  single  crystals  or  borosilicate  glass.  It  is  demonstrated  that  all  relevant  optical 
and  thermal  properties  of  the  dosimeters  can  be  determined  by  these  methods  and  that,  based  on 
this  knowledge,  the  expected  thermoluminescence  response  of  a  given  configuration  can  be 
simulated  as  a  function  of  a  number  of  experimental  parameters. 


I.  INTRODUCTION 

In  recent  years  the  feasibility  of  heating  thermolu¬ 
minescent  materials  with  infrared  laser  beams  has  been  dem¬ 
onstrated.  This  new  thermal  stimulation  method  has  re¬ 
ceived  particular  attention  in  thermoluminescence 
dosimetry  (TLD)  of  ionizing  radiation  because  it  holds  pro¬ 
mise  as  a  solution  to  a  number  of  problems  associated  with 
the  measurement  of  small  doses  of  nonpenetrating  radiation 
such  as  low-energy  fi  rays  as  well  as  knock-on  protons  pro¬ 
duced  by  fast  neutrons  in  hydrogenous  radiators.5-7  Even 
the  imaging  of  spatial  dose  distributions  of  y  rays  and  x  rays 
has  been  contemplated 1,2  and  experimentally  demonstrat¬ 
ed.8 

The  successful  development  of  practical  laser-heated 
TLD  readers  and  thermoluminescence  imaging  devices,  re¬ 
quires  detailed  knowledge  of  the  heat  transfer  mechanism 
from  the  laser  beam  to  the  thermoluminescent  sample  and 
the  resulting  spatiotemporal  temperature  distribution.  Then 
one  may  compute  and  predict  the  thermoluminescence  re¬ 
sponse  (glow  curve)  of  the  sample  dosimeter  for  a  given  laser 
heating  beam.9  The  effects  of  such  interrelated  design  pa¬ 
rameters  as  laser  power  density,  time  of  laser  exposure,  laser 
beam  diameter  and  intensity  profile,  sample  configuration 
and  composition,  and  the  expected  incandescent  back¬ 
ground  emission  can  be  computer-simulated  and  the  result¬ 
ing  glow  curve  compared  with  experimental  results. 

The  list  of  experimental  parameters  is  narrowed  down 
to  the  few  that  presently  appear  important  for  practical  ap¬ 
plications.  For  example,  the  laser  power  of  cost-effective 
modem  rf  or  dc  discharge  excited  C02  waveguide  lasers  or 
rf-excited  unstable  resonator  lasers  is  limited  to  less  than  100 
W.  Laser  beam  profiles  are  normally  Gaussian  intensity  dis¬ 
tributions  or  the  characteristic  annular  cross  section 
(“halo”)  that  is  possible  with  an  unstable  resonator  cavity. 
Uniform  beam  profiles  are  highly  desirable,  but  all  known 
methods  to  produce  such  beams  appear  to  be  impractical  in 


TLD  applications  because  of  the  associated  power  loss  or  the 
required  custom  fabrication  of  expensive  special  optics. 10 

Suitable  dosimeter  configurations  are  thin  layers  of  the 
TLD  phosphor,  either  in  the  form  of  a  self-supporting  film 
composed  of  a  high  temperature  polymer  mixture  with  phos¬ 
phor  particles,  or  applied  to  a  thin  substrate.  Both  may  take 
the  form  of  discrete  circular  spots  whose  diameter  is  smaller 
than  the  full  width  at  1/e  power  of  the  laser  beam  or  may  be 
prepared  as  continuous  layers  whose  area  is  much  larger 
than  tne  laser  beam  diameter.  However,  even  the  widely 
used  square  TLD  chips  (e.g.,  the  configurations  supplied  by 
Harshaw/Filtrol")  or  round  pellets12  can  be  successfully 
heated  with  a  CO:-laser  beam.13 

Small  discrete  dosimeters  can  be  heated  quite  uniformly 
with  Gaussian  beams  provided  the  laser  photons  are  not  too 
strongly  absorbed  (uniform  heating  in  the  direction  of  laser 
propagation).  Strongly  absorbing  dosimeter  layers  (surface 
heating)  must  be  thin  enough  that  the  thermal  response 
time14  of  the  sample  is  much  shorter  than  the  total  heating 
time  required  to  release  all  trapped  carriers.  A  necessary 
condition  for  a  spatially  uniform  temperature  rise  in  both 
cases  is,  of  course,  that  the  diameter  of  the  discrete  spot  is 
much  smaller  than  the  diameter  of  the  laser  beam.  Obvious¬ 
ly,  a  considerable  fraction  of  the  laser  beam  energy  is  wasted 
in  this  mode  of  operation  which,  for  this  reason,  will  prob¬ 
ably  be  considered  only  in  special  situations. 

Continuous  dosimeter  layers  are  of  practical  interest  be¬ 
cause  of  the  ease  with  which  they  can  be  fabricated  and  be¬ 
cause  of  their  potential  usefulness  in  a  number  of  the  new 
applications  mentioned  above.  A  localized  Gaussian  beam 
can  be  employed  to  rapidly  heat  a  small  spot  on  the  TLD 
layer,  yielding  a  characteristic  glow  curve  that  is  the  result  of 
the  nonuniform  spatiotemporal  temperature  distribution 
T\x,yj,t )  produced  in  this  way.  This  type  of  thermolumines¬ 
cence  response  is  completely  different  from  the  conventional 
glow  curves  obtained  by  uniform  contact  heating  (an  exam¬ 
ple  is  shown  in  Fig.  9  below).  It  is  nevertheless  very  useful  in 
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dosimetric  and  imaging  applications  of  thermolumines¬ 
cence.15 

In  this  paper  we  present  the  theory  of  the  thermolu- 
minescence  response  curves  generated  with  Gaussian  laser 
beam  profiles  in  a  number  of  dosimeter  configurations  that 
are  of  practical  importance  in  thermoluminescence  dosi¬ 
metry.  The  theoretical  approach  is  verified  by  comparison 
with  experimental  results.  In  addition,  experimental  meth¬ 
ods  are  presented  for  the  determination  of  all  relevant  ther¬ 
mal  and  optical  properties  of  the  dosimeter  materials.  Thus, 
the  thermoluminescence  response  of  a  given  type  of  dosi¬ 
meter  is  completely  characterized  and  the  effects  of  such 
design  parameters  as  layer  thickness,  laser  power  and  beam 
size,  preannealing  temperature  and  duration,  substrate  ma¬ 
terial  and  its  thickness  can  all  be  assessed  by  computer  simu¬ 
lation.  The  special  case  of  semi-infinite  LiF  slab  has  been 
reported  previously.9 

II.  THERMOLUMINESCENCE  KINETICS  FOR 
NONUNIFORM  SPATIOTEMPORAL  TEMPERATURE 
DISTRIBUTIONS 

A.  General  remarks 

Calculations  of  thermoluminescence  glow  curves  on 
the  basis  of  simple  electron  kinetic  models  and  for  spatially 
uniform  heating  pose  no  principle  problems.  While  only  spe¬ 
cial  cases  yield  analytical  solutions  for  the  thermolumines¬ 
cence  intensity  as  a  function  of  time  or  temperature,  numeri¬ 
cal  solutions  can  readily  be  generated  for  the  coupled 
nonlinear  rate  equations  which  describe  a  given  trap  mod¬ 
el.  1617  The  extreme  heating  rates  possible  with  lasers  have 
been  shown  not  to  invalidate  the  electron-statistical  founda¬ 
tion  (Shockley-Read  statistics1819)  of  these  phenomenologi¬ 
cal  theories.4 

Usually  a  number  of  different  physically  plausible  elec¬ 
tron-kinetic  models  will  give  satisfactory  fits  of  experimental 
glow  curves.  This  general  lack  of  uniqueness  renders  trap- 
level  spectroscopy  by  monitoring  thermally  simulated  elec¬ 
tron-kinetic  relaxation  phenomena  a  nontrivial  task.1617 
Spatially  nonuniform  heating  significantly  increases  the  re¬ 
quired  computational  effort  and  further  complicates  this  sit¬ 
uation.  For  this  reason,  laser  heating  techniques  are  ex¬ 
plored  herein  for  the  sole  purpose  of  assessing  their  utility  in 
the  dosimetry  of  ionizing  radiation.  No  attempt  is  made  to 
verify  a  given  simple  trap  model  or  to  extract  new  methods 
which  might  facilitate  trap-level  spectroscopy.  Thus,  all  that 
is  required  for  our  purposes  is  one  reasonably  plausible  mod¬ 
el  description  for  the  thermoluminescence  kinetics  of  a  suit¬ 
able  phosphor  heated  with  a  laser. 

The  calculations  presented  here  are  based  on  simple 
Randall- Wilkins  first-order  kinetics. 16  As  the  example  for  a 
typical  phosphor,  we  choose  the  most  widely  used  dosimetry 
material  LiF:Mg,Ti  (TLD- 100). 12  Experimental  fits  of  its 
prominent  glow  peaks  (costumarily  labelled  peaks  2-5)  have 
been  reported  by  a  number  of  authors,  most  recently  by 
McKeever20  and  by  Vana  and  Ritzinger.21  The  fact  that  the 
numerical  fitting  parameters,  notably  the  trap  depths  and 
the  frequency  factors,  differ  significantly  in  these  papers 
only  underlines  the  problems  associated  with  curve  fitting  as 
a  method  to  obtain  basic  physical  information  from  thermo¬ 


luminescence  experiments  alone.  We  have  arbitrarily  select¬ 
ed  McKeever’s  data  for  our  investigations. 

B.  Computational  approach 

Having  selected  a  group  of  glow  peaks  in  TLD- 100  and 
a  suitable,  albeit,  not  unique,  fit  for  modeling  these  p~»v» 
under  various  experimental  conditions,  the  main  task  is  the 
calculation  of  the  time  evolution  of  the  temperature  (“heat¬ 
ing  program”)  for  each  volume  element  dxdydz  at  location 
[x,yj:)  of  the  sample. 

A  rich  literature  exists  in  the  field  of  laser  heating  of 
materials.  For  example,  laser  annealing  of  semiconductors 
has  become  an  important  manufacturing  tool  in  the  elec¬ 
tronics  industry.22'27  Laser  absorption  calorimetry  of  very 
transparent  optical  materials  also  requires  the  knowledge  of 
the  time-dependent  temperature  distributions  produced  in 
the  sample.28  Important  insight  in  the  laser  damage  proper¬ 
ties  of  high-power  metallic  and  dielectric  mirrors  and  other 
optical  components  has  been  gained  by  studying  the  theory 
of  laser  heating.29  30  Comprehensive  treatments  of  heat  con¬ 
duction  in  solids  and  associated  boundary  value  problems 
are  the  texts  by  Carslaw  and  Jaeger31  and  Ozisik.32  However, 
despite  an  extensive  bibliography  on  laser  heating,  the  prob¬ 
lem  of  heating  thermoiuminescence  dosimeters  with  a  Gaus¬ 
sian  beam  turned  out  to  be  a  unique  new  case.  It  is  perhaps 
most  closely  related  to  the  one  solved  by  Bernal2*  who  con¬ 
sidered  Gaussian  beam  heating  of  extremely  weakly  absorb¬ 
ing  cylinders  and  semi-infinite  slabs.  The  laser  beam  is,  in 
first  approximation,  unattenuated  when  traversing  the  en¬ 
tire  thickness  of  the  slab.  Unfortunately,  this  special  case  is 
of  little  interest  here.  Weakly  absorbing  thermoiumines¬ 
cence  dosimeter  layers  are  not  suited  for  practical  applica¬ 
tions  because  of  inefficient  use  of  costly  beam  energy.  Strong 
to  moderate  beam  attenuation  in  the  active  layer  is  charac¬ 
teristic  for  laser  heating  in  thermoiuminescence  dosimetry. 
However,  this  problem  can  be  solved  analytically  by  the 
Green’s  function  method  under  the  assumption  that  the  lay¬ 
er  is  semi-infinite  in  radial  direction  (local  heating  of  a  small 
spot  on  a  large  area)  and  that  no  significant  heat  loss  occurs 
due  to  convection  and  radiation.9  Convection  can  easily  be 
eliminated  by  evacuation  of  the  sample  chamber  and  radia¬ 
tive  losses  are  indeed  small  considering  the  fact  that  the  max¬ 
imum  temperatures  of  interest  do  not  exceed  700  K  and  the 
heating  times  are  limited  to  a  few  hundred  milliseconds.  For 
these  assumptions  to  be  valid,  the  rate  of  energy  absorption 
from  the  laser  beam  must  be  significantly  greater  than  that  of 
radiative  energy  loss. 

The  resulting  spatiotemporal  temperature  distribution 
is  cylindrically  symmetric  around  the  beam  axis.  To  calcu¬ 
late  the  total  thermoluminescent  emission  as  a  function  of 
time  after  onset  of  the  laser  exposure,  the  time-dependent 
contributions  from  volume  elements  "Lnr&rAi  are  summed 
by  numerical  integration  over  the  radial  coordinate  r  and  the 
direction  of  laser  beam  propagation  z.  Because  these  heating 
rates  are  nonlinear  in  time,  analytic  solutions  of  the  electron 
kinetic  rate  equations  for  each  of  these  infinitesimal  volume 
elements  and  their  individual  unique  heating  rate  are  un¬ 
available  even  for  simple  first-order  kinetics.  Therefore,  the 
rate  equations  must  be  solved  numerically  on  the  basis  of  the 
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individual  temperature  evolution  of  each  of  these  volume 
elements  and  their  contribution  to  the  emission  intensity  I  (r ) 
must  be  summed  up. 

The  typical  situation  encountered  in  heating  a  small 
area  of  a  continuous  semi-infinite  layer  of  thickness  L  by  a 
Gaussian  laser  beam  of  intenstiy 

/  (rj)  =  I0  expi  -  r2/ w 2  -  fiz)  ( 1 ) 

is  schematically  depicted  in  Fig.  1  The  symbols  chosen  are 
the  full  width  2w  at  i/e  peak  power  of  the  beam,  the  absorp¬ 
tion  coefficient  fi,  the  radial  distance  r  from  the  axis,  and  the 
distance  z  from  the  front  surface  upon  which  the  beam  is 
incident.  If  the  beam  is  turned  on  at  t  =  0,  the  temperature 
distribution  T\rj,t )  evolves  as  a  consequence  of  energy  ab¬ 
sorption  and  thermal  diffusion,  resulting  in  a  spatially  re¬ 
solved  time-dependent  thermoluminescence  emission  pat¬ 
tern  ITL\r,L,t )  similar  to  the  special  case  shown  in  Fig.  2. 
These  photographs  were  obtained  by  heating  a  highly  ab¬ 
sorbing  1-mm-thick  glass  slab  whose  back  side  at  z  =  L  was 
coated  with  a  50-^m-thick  layer  of  ZnS:Cu  powder.  Hefe  the 
thermoluminescence  emission  ITL(rj,t )  is  not  due  to  direct 
laser  heating  of  the  phosphor  itself.  Instead,  the  thermolu¬ 
minescent  layer  is  heated  by  temperature  diffusion  through 
the  glass.  Similar  experiments  were  described  previously.2 
After  appropriate  calibration  such  a  thermoluminescent  lay¬ 
er  can  be  used  to  measure  the  time  evolution  of  isotherms 
during  laser  heating  of  materials.  For  example,  the  diameter 
of  the  dark  center  spot  in  Figs.  2(b)  and  2(c)  corresponds  to 
the  438  K  isotherm  because  the  thermoluminescence  emis¬ 
sion  of  ZnS:Cu,  measured  by  conventional  uniform  slow 
heating  at  a  constant  heating  rate  of  15.3  K/s,  ceases  at  this 
temperature.  Note  that  the  average  heating  rate  in  the  laser 
experiment  (Fig.  2)  is  about  the  same  as  can  be  discerned 
from  the  appearance  of  the  dark  center  spot  at  about  10  s. 
Other  isotherms  may  be  selected  as  well.  For  example  the 
peak  emission  (region  of  maximum  brightness  in  Fig.  2)  cor¬ 
responds  to  approximately  .'31  K.  Still  other  emission  inten¬ 
sities,  measured  relative  to  the  peak  emission,  may  be  cho- 


FIG.  1.  Schematic  representation  of  a  Gaussian  laser  beam  profile  exposing 
the  surface  iz  =  0)  of  a  slab  of  thickness  L 


FIG.  2.  Thermoluminescent  emission  pattern  obtained  from  a  50-/im-thick 
uniform  ZnS  Cu  layer  on  a  1-mm-thick  glass  slab  |?5  x  25  mm*).  The  sam¬ 
ple  was  heated  by  a  cw  C02  laser  beam  of  4  3  W  and  5.4-mm  full  width  at  1  / 
e  peak  power.  The  beam  was  directed  perpendicularly  at  the  center  of  the 
glass  slab  on  the  side  opposite  to  the  phosphor  layer.  The  photographs  were 
taken  at  the  indicate  j  times  after  onset  of  the  laser  exposure. 


sen.  Actually,  the  entire  temperature  distribution  present  at 
a  given  time  after  onset  of  heating  can  be  determined  by 
comparison  of  the  normalized  emission  intensities  in  Fig.  2 
with  those  of  the  conventional  glow  curve  measured  at  ap¬ 
proximately  the  same  heating  rate. 

In  dosimetry  applications  of  laser  heating  the  emission, 
integrated  over  the  total  sample  volume,  is  measured  with  a 
photomultiplier  tube  without  spatial  resolution: 

A-l(')  =  j jdrdzl„(rj,t).  (2) 

Here  Iv  is  the  thermoluminescence  emission  intensity  per 
unit  volume. 

As  a  first  step  in  the  theory  of  this  TL  emission  we  calcu¬ 
late  T\rj,t )  by  solving  the  thermal  diffusion  equation9 

V27>,z,r)  +  k  ~'g(rj,t)  =  a~'dT/dt  (3) 

for  the  following  boundary  conditions. 

(a)  The  slab  is  infinite  in  radial  direction  and  has  a  thick¬ 
ness  L. 

(b)  No  heat  loss  occurs  at  the  boundaries  z  =  0  and 
z  =  L,  i.e. ,  dT /dz  =  0.  This  assumption  is  valid  in  vacuum  in 
the  absence  of  radiative  heat  loss. 

The  initial  condition  is  T[rj,t)=  T0,  i.e.,  room  tem¬ 
perature  when  the  laser  beam  is  turned  on  at  t  =  0.  Then  the 
temperature  obtained  from  Eq.  (3)  is  the  actual  temperature 
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increase  above  T0.  The  notation  used  in  Eq.  (3)  is  as  follows: 
k  =  thermal  conductivity,  a  =  k  /pc  =  thermal  diffusivity, 
p  =  mass  density,  c  =  specific  heat,  girj.t )  =  laser  power 
per  unit  volume  present  at  position  (r,z)  during  laser  expo¬ 
sure  (source  function). 

The  case  of  temperature-dependent  k  and  a  is  discussed 
in  Appendix  A. 

Correcting  for  reflective  loss  on  the  front  surface  at 
z  =  0,  the  source  function  g  becomes 

=  ui  1  —  R  )/0  exp|  —  r/w 2  —  az\  for  :>0 

fflW)  n  f  ,  n  i4> 

=  0  for  t  <  0, 

where  R  is  the  reflectivity. 

We  solve  Eq.  (3)  with  the  initial  and  boundary  condi¬ 
tions  stated  above  by  the  Green’s  function  technique.32  The 
general  solution  is 

AT  =T [rj,t )  -  T0  =  f  drf  r'dr'  f  dz' 
k  Jr  =  o  Jr1  =  0  Jl'  =•  0 

=  [G{rj,r,r’j'r)g{r'j'r)],  (5) 

where  G  is  the  following  Green’s  function: 

G  [rj,t;r' J ,t\ 


-f' 

Je 


d0e 


-  aB:\t 


WofirVotfr') 


X{j[Jr\'Le  aVm['  n  coslnmz)  cos(i7mr')j.  (6) 

Here/0  is  the  zero-order  Bessel  function  and  rjm  =  mirL  ~ 
After  integration  over  r’£,  and  /?,  one  obtains  with 
P  —  f q  I  (r,0)2irrdr  —  mvzI0  (total  power  incident  on  the 
front  surface  of  the  thermoluminescent  layer  at  z  =  0)  the 
temperature  distribution9: 

at=[- - dr- - - 

kLir  Jr-o  4a(f  —  r)  +  w 


X 


30 

+  2I 

m  =  t 


e  -  avU<  -  r 

1  +iv„/p)2 


X 1 1  -  e  -  uL 


cos  mn) 


(7a) 


Note  again  that  Eq.  (7a)  is  the  temperature  in  excess  of  the 
initial  temperature  T0.  The  properties  of  AT  (rj,t )  in  Eq.  (7a) 
can  be  readily  appreciated  by  discussing  two  limiting  cases 
for  the  absorption  coefficient  p.  Weak  absorption,  that  is 
pL<  1,  means  uniform  energy  absorption  in  the  z  direction. 
Thez-dependent  term  in  the  second  bracket  vanishes  and  the 
remaining  factor  (1  —  e  ~  )  is  simply  the  fraction  of  the 

energy  absorbed  from  the  transmitted  Gaussian  beam.  Dif¬ 
fusion  broadens  the  beam  profile  with  increasing  time.  Very 
strong  absorption,  that  is  pL>  1,  means  only  a  thin  layer  on 
the  surface  at  z  =  0  is  directly  heated  and  bulk  heating  oc¬ 
curs  only  by  thermal  diffusion.  This  case  is  of  interest  for 
certain  practical  dosimeter  configurations. 

The  actual  computer  simulations  make  use  of  the  heat¬ 
ing  rate  of  each  volume  element  lirrArAz,  calculated  from 
Eq.  (7a): 


( l  _  p  \np  e  -  +  •*) 

dT/dt  =  - - — 

kLir  4 at  4-  w 


X 


(Vm/Pf 


X  ( 1  —  e  “L  cos  mir)  cosir) mz)J.  (7b! 

The  four  first-order  differential  equations  which  de¬ 
scribe  the  independent  processes  of  peaks  2-520  were  solved 
together  with  the  heating  rate  of  Eq.  (7b)  by  the  Runge- 
Kutta  method  for  each  cell  ArAz.  Sufficient  precision  was 
obtained  with  A  r  =  2w/50,  Az  =  L  / 50  and  A  t  =  1 0  ~ 3  x  the 
total  exposure  time.  The  sum  of  the  relative  peak  heights211 
was  used  as  a  measure  of  the  total  trapped  electron  popula 
tion  and  each  peak  height  was  normalized  by  this  sum  to 
provide  the  initial  conditions  for  the  TL  differential  equa¬ 
tions.  Annealing  of  the  lowest  temperature  peak  no.  2  or  any 
other  is  easily  simulated  by  putting  its  relative  height  to  zero. 


C.  Experimental  procedure 

The  facility  shown  schematically  in  Fig.  3  was  con¬ 
structed  for  the  measurement  of  laser  stimulated  thermolu¬ 
minescence  response  curves. 

A  radio-frequency  excited  C02  laser  (Laakmann  Elec¬ 
tro-Optics,  rated  at  4  W)  delivers  a  cw  beam  of  1.3-mm  full 
width  at  1/e  maximum.  Together  with  all  the  required  opti¬ 
cal  elements  and  positioners,  this  laser  is  mounted  on  a  stur¬ 
dy  aluminum  table. 

The  electromechanical  shutter  is  the  first  optical  ele¬ 
ment  in  the  laser  beam  path.  To  avoid  thermal  distortion  of 
the  shutter  blades,  they  are  gold  plated  for  high  reflectivity. 
During  the  closed  position  of  the  shutter,  the  laser  beam  is 
reflected  from  the  blades  into  a  beam  dump  which  simply 
absorbs  and  dissipates  the  laser  power.  The  shutter  has  a  dc 
power  supply  which  is  housed  inside  the  electronic  control 
enclosure.  The  electronic  controlled  shutter  is  open  only  to 
heat  a  sample  on  the  sample  holder  or  to  measure  by  a 
powermeter  the  laser  power  delivered  to  the  sample.  The 
powermeter  is  mounted  on  a  shaft  which  slides  through  the 
light  tight  box. 

Following  the  laser  beam  path  beyond  the  shutter,  a  5  x 
beam  expander  is  mounted  on  the  reader  table.  It  can  be 
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FIG.  3.  Schematic  representation  of  the  laser-heated  thennolummescence 
reader. 
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adjusted  to  produce  various  beam  diameters  (1.5-6  mm  full 
width)  entering  the  focusing  cube.  The  focusing  cube  deflects 
the  beam  by  90*  through  an  AR-coated  germanium  window 
and  causes  it  to  converge  on  the  TL  sample  in  the  light  tight 
box.  The  germanium  window  is  transparent  at  the  10.6-/xm 
laser  wavelength,  but  opaque  to  visible  light.  The  beam  posi¬ 
tion  on  the  TL  sample  is  determined  by  moving  the  sample 
holder  with  two  connected  x-y  translation  stages.  The  trans¬ 
lation  stages  have  two  actuators  which  are  driven  by  the 
Newport  Corporation  model  350M  motor  drive  system  and 
give  an  accuracy  of  the  beam  position  up  to  0. 1  /xm. 

The  TL  sample  is  heated  as  it  absorbs  this  infrared  radi¬ 
ation  from  the  laser.  The  high-power  density  of  the  focused 
laser  beam  causes  the  selected  spot  on  the  TL  sample  to  be 
heated  rapidly,  while  the  remainder  of  the  sample  remains 
cool. 

The  TL  emitted  by  the  sample  passes  the  sample  sub¬ 
strate  and  reaches  into  the  photomultiplier  tube  (PMT)  be¬ 
low.  The  PMT  is  protected  from  two  types  of  accidental 
exposure,  one  due  to  the  laser  by  a  glass  window  which  will 
absorb  the  10.6-/xm  laser  beam,  another  one  due  to  the  visi¬ 
ble  light  I  when  the  sample  chamber  is  open)  by  a  thick  alumi¬ 
num  slide.  The  PMT  is  an  EMI  type  9924B  with  a  l-in.-diam 
bialkali  photocathode.  The  current  signal  from  the  PMT  is 
digitized,  stored,  and  displayed  by  a  Data  Precision  6000 
Universal  Waveform  Analyzer  (transient  digitizer). 

The  laser  TLD  reader  was  designed  to  permit  examina¬ 
tion  of  a  number  of  experimental  parameters  related  to  laser 
heating  TLD,  e  g.,  the  laser  power  density  that  the  sample 
receives,  and  the  beam  diameter  which  is  set  by  mechanically 
adjusting  the  beam  expander.  Laser  power  and  heating  time 
are  controlled  by  settings  in  the  electronic  control  unit  (Fig. 
4). 

The  timing  circuits  in  the  electronic  control  unit  consist 


FIG  4  Block  diagram  of  the  electronic  control  system  for  the  apparatus 
shown  in  Fig.  3 


of  four  LM555  timers  that  generate  four  timing  pulses.  The 
first  pulse  drives  the  cycle  indicator  and  limits  read  cycles  to 
one  per  second.  The  second  pulse  is  used  to  turn  the  laser  off 
and  open  the  shutter.  The  third  pulse  turns  the  laser  on  for  a 
precise  preselected  timing  interval.  The  transient  digitizer  is 
triggered  at  the  beginning  of  this  pulse,  and  the  laser  is 
turned  off  at  its  end.  The  fourth  pulse  triggers  the  closing 
mechansim  of  the  shutter.  The  laser  power  during  the  third 
timing  pulse  is  determined  by  modulation  control.  When¬ 
ever  all  four  of  these  pulses  are  off,  the  laser  operates  in  a  cw 
mode  at  maximum  output  power  to  maintain  a  stable  operat¬ 
ing  temperature. 

The  modulation  circuit  consists  of  a  10-MHz  oscillator 
that  triggers  a  variable  length  pulse  generator.  The  modula¬ 
tion  control  is  used  to  vary  the  pulse  width  of  the  resulting 
pulse  train,  permitting  the  laser  output  power  to  be  varied 
over  a  range  from  zero  to  full  power.  The  modulation  driver 
amplifies  the  modulation  signal  to  the  6-V  pulse  height  re¬ 
quired  for  the  laser  rf  power  supply.  The  shutter  control 
signal  and  the  data  signal  are  both  logic  level  signals.  The 
74LS  series  TTL  logic  circuits  are  employed  for  all  control 
logic. 

As  will  be  demonstrated  below  (see,  for  example,  Fig. 
12),  the  shape  of  the  thermoluminescence  response  curve  de¬ 
pends  critically  on  the  laser  beam  profile  in  the  sample  plane 
at  z  =  0.  For  this  reason  we  have  given  its  precise  measure¬ 
ment  particular  attention,  using  the  so-called  scanning  knife 
edge  method.33 

A  photodetector  measures  the  laser  power  not  occluded 
by  a  knife  edge,  scanning  through  the  beam,  as  a  function  of 
position.  The  spatial  profile  is  then  given  by  the  first  deriva¬ 
tive  of  the  photedetector  signal  with  respect  to  the  scan  di¬ 
rection.  A  necessary  requirement  is  that  the  mathematical 
form  of  the  profile  be  separable  in  the  coordinates  perpendic¬ 
ular  to  and  parallel  to  the  scan  directions  as  in  the  case  of  a 
Gaussian. 

To  check  for  deviations  from  circular  symmetry,  we 
have  always  scanned  in  two  perpendicular  directions.  For 
cw  lasers  this  procedure  is  readily  executed  by  scanning  at  a 
constant  rate  and  analog  differentiation  with  an  appropriate 
RC  circuit.34  However,  for  infrared  lasers  such  as  the  one 
used  in  our  work,  thermal  detectors  are  commonly  em¬ 
ployed  to  monitor  the  beam  power.  Because  of  their  satura¬ 
tion  characteristics  they  are  generally  not  suitable  for  scan¬ 
ning  knife-edge  applications.  Modifications  of  this  method 
are  required.  They  are  shown  schematically  in  Fig.  5.  The  cw 
beam  is  chopped  with  the  electromechanical  shutter3  driven 
by  an  oscillator.  The  oscillator  output  is  used  simultaneously 
as  a  clock  for  a  transient  digitizer  (Data  Precision  model  DP 
6000)  to  synchronize  the  sampling  interval  with  the  occur- 
rance  of  the  signal  pulse  from  a  piezoelectric  detector 
( Barnes  Engineering  model  3  50  PZT).  In  our  case  the  shutter 
is  opened  ever  73.5  ms  for  3.5  ms.  The  recorded  detector 
signal  during  a  complete  scan  has  the  appearance  of  a  digital¬ 
ly  recorded  cw  signal,  exhibiting  only  the  slow  temproal 
structure  due  to  increasing  occlusion  of  the  beam  by  the 
knife  edge.  This  signal  is  then  differentiated  mathematically 
by  the  microprocessor  of  the  DP  6000. 

Examples  of  beam  profiles  obtained  in  this  fashion  are 
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KNIFE-  PZT 

SHUTTER  LENS  EDGE  OETECTOR 


FIG.  5  Schematic  arrangement  for  knife-edge  scanning  of  the  CO,  laser 
beam  profile.  Note  that  the  position  of  the  knife-edge  in  the  sample  plane  is 
shown  here  vertically  because  the  focusing  cube  isee  Fig.  3)  is  omitted  for 
simplicity. 


shown  in  Fig.  6.  The  combinaton  of  the  beam  expander  and 
5-in.  focal  length  focusing  cube  afforded  variations  of  the 
laser  spot  size  in  the  sample  plane  between  0.3  and  2.22  mm. 
The  long  focal  length  of  the  cube  and  the  maximum  laser 
beam  penetration  depth  in  our  samples  of  less  than  300  /itn 
implies  that  the  laser  beam  profile  is  independent  of  the  z 
direction. 

In  the  following  sections  we  present  experimental  re¬ 
sults  obtained  with  a  number  of  different  dosimeter  configu¬ 
rations  and  comparisons  with  theoretical  calculations. 

III.  RESULTS  AND  DISCUSSION 

Case  1:  LiF.  Ti.  Mg  chip 

Laser  heating  of  microcrystalline  hot-pressed  TLD-100 
dosimeter  material  produced  by  the  Harshaw/Filtrol  Cor¬ 
poration11  is  of  interest  as  it  is  widely  employed  in  routine 
personnel  dosimetry.  From  a  theoretical  point  of  view  this 
case  is  important  because  of  the  rather  modest  absorption 
coefficient  ( /*  =  40  cm" ')  of  LiF  for  10.6 ptm  photons.  The 
laser  beam  penetrates  250  n m  (1/e  attenuation  depth)  into 
the  900-/m-thick  chip.  As  it  turns  out,  none  of  the  previous¬ 
ly  reported  solutions  of  the  heat  flow  equation  [Eq.  (3)]  cov¬ 
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FIG  6.  Example  of  a  typical  nearly  Gaussian  spatial  profile  obtained  for  the 
C03  laser  beam  at  r  =  0  (see  Fig  l|. 


ers  this  situation  in  conjunction  with  a  Gaussian  heating 
beam. 

Unfortunately  only  small  chips  (3 X 3x0.9  mm3)  are 
commercially  available  and,  therefore,  care  has  to  be  taken 
not  to  violate  the  assumption  of  a  semi-infinite  medium  (see 
Sec.  II)  when  comparing  computed  with  experimental  re¬ 
sults.  In  practice  this  can  be  done  by  choosing  a  beam  diame¬ 
ter  that  is  much  smaller  than  the  width  of  the  sample  and  by 
limiting  the  total  exposure  of  the  laser  heating  beam  to  a  time 
during  which  the  sample  edges  have  not  yet  been  heated  by 
thermal  diffusion.  In  order  to  reach,  during  this  time,  in  the 
center  of  the  exposed  area  the  highest  temperature  peak 
(peak  no.  5),  the  power  density  in  the  beam  center  must  ex¬ 
ceed  a  minimum  value.  In  essence  this  requirement  means 
that  a  substantial  part  of  the  chip  volume  must  be  heated  to 
above  the  maximum  temperature  for  which  TL  emission  will 
still  occur  before  the  edge  areas  experience  any  temperature 
increase  by  thermal  diffusion.  We  have  performed  these  ex¬ 
periments  successfully  with  beam  diameters  around  0.084 
cm  and  approximately  5-W  power.  Under  these  conditions  a 
300-ms  exposure  time  raised  the  temperature  at  the  edges  by 
only  a  few  degrees  K  above  room  temperature.  Continued 
exposure  to  the  heating  beam  beyond  this  time  resulted  in 
heat  pileup  in  the  chip  and  eventual  increase  of  the  tempera¬ 
ture  of  the  total  sample  to  above  700  K,  a  situation  not  cov¬ 
ered  by  the  theory  presented  in  Sec.  II  B.  However,  these 
experiments  are  perhaps  important  as  they  indicate  that  en¬ 
tire  TLD-100  or  similar  chips  can  indeed  be  heated  with  an 
appropriate  laser  beam  so  as  to  measure  the  total  thermolu- 
minescence  emission.  In  fact,  we  have  been  able  to  heat  these 
chips  to  600  K  within  about  2  s  with  a  2.5-mm-diam  beam  of 
approximately  3-W  power. 

The  temperature  distribution  established  at  the  end  of 
the  exposure  1 300  ms)  of  a  LiF  chip  to  a  0.084-cm  full  width 
at  l/e  maximum  of  4.93  W  power  is  depicted  in  Fig.  7.  Iso¬ 
therms,  calculated  from  Eq.  (7a),  are  shown  for  temperatures 
corresponding  to  the  peak  temperatures  of  peaks  2-5  (383, 
421,  457,  and  483  K,  respectively).  The  volume  inside  the 
483  K  isotherm  is  heated  to  temperatures  for  which  thermo¬ 
luminescent  emission  from  peaks  2-4  has  ceased  at  300  ms. 
Figure  7  illustrates  the  nonuniformity  of  both  the  tempera¬ 
ture  distribution  and,  albeit  indirectly,  the  thermolumines¬ 
cent  emission  pattern. 

We  have  made  no  attempts  in  our  experiments  to  mea¬ 
sure  the  spatial  distribution  of  the  time-dependent  thermolu¬ 
minescence  response  of  the  TLD-100  chip  to  the  Gaussian 
heating  beam.  Instead  the  total  emission  /XL(r )  was  moni¬ 
tored  with  a  photomultiplier  tube  and  compared  with  com¬ 
putations  performed  with  the  theory  outlined  in  Sec.  II  B. 
The  results  are  presented  in  Fig.  8  for  a  TLD-100  chip  that 
was  preannealed  to  373  K  for  10  min  in  order  to  remove  peak 
2.  This  type  of  thermal  treatment  after  exposure  to  the  ioniz¬ 
ing  radiation  and  before  measuring  the  thermoluminescence 
glow  curve  is  customary  in  dosimetry  applications  of  this 
material.  The  agreement  between  the  calculated  and  experi¬ 
mental  glow  curve  is  rather  satisfying  in  light  of  the  uncer¬ 
tainties  concerning  the  exact  electron-kinetic  model  descrip¬ 
tion  of  the  thermoluminescence  emission  from  this  material 
and  the  fact  that  the  TLD-100  chip  is  not  a  single  crystal  as 
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FIG.  7.  Temperature  distribution  established 
after  heating  a  small  LiF  slab  of  area  3x3  mm2 
and  thickness  L  =  0.9  mm  for  300  ms.  A  Gaus¬ 
sian  laser  beam  of  4.93-W  power  and  0.84-mm 
full  width  impinges  on  the  center  of  the  slab  at 
z/L  =  0.  The  calculations  were  performed  with 
Eq.  I7b|  and  k  =  0.0275  W/cm  K,  c  =  1  56 
W  s/g  K.  p  =  2.633  g/cm\  and  an  absorption 
coefficient  for  10/rm  of/i  =  40  cm  ~ '. 
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assumed  in  the  calculations.  An  additional  uncertainty 
stems  from  the  rather  strong  temperature  dependence  of  the 
thermal  conductivity  k.  According  to  Men  et  al.,iS  k  de¬ 
creases  with  increasing  temperature  approximately  as  T  ~ '. 
The  calculated  curve  shown  in  Fig.  8  was  obtained  with 
k  =  0.0275  W/cm  K,  which  is  appproximately  one  third  of 
the  value  measured  for  a  single  crystal  at  room  temperature. 
We  have  also  performed  the  calculation  with  an  exact  solu¬ 
tion  of  the  heat  flow  equation  with  k[T)  —  kQTJT  using 
Kirchhoff  transformation  techniques  (see  Appendix  A). 
However,  no  significant  improvement  in  the  agreement 
between  measured  and  calculated  thermoluminescence  re¬ 
sponse  curves  was  obtained,  perhaps  indicating  that  unac¬ 
counted  for  experimental  factors  influenced  the  results  more 
than  an  improved  theory  of  heat  conduction.  We  suspect 
that  the  microcrystalline  “ceramic”  nature  of  the  sample 
and  the  roughness  of  the  unpolished  surface  caused  some 
deviations  from  the  implicit  assumption  that  the  optical  and 
thermal  properties  of  this  material  are  very  close  to  those  of 
single-crystalline  LiF.  Scattering  of  the  laser  light  probably 
results  in  some  broadening  of  the  Gaussian  beam  profile  in¬ 
side  the  chip,  and  the  thermal  conductivity  may  be  lower 
than  that  of  a  single  crystal.  In  addition,  it  is  well  known 
that,  in  order  to  reproduce  the  exact  relative  peak  heights  of 
peaks  2-5  of  the  TLD-100  phosphor  as  used  in  Ref.  20,  the 
annealing  procedure  must  be  reproduced  exactly.12  Slight 
deviations  might  have  been  present  as  compared  to  those 
used  by  McKeever,20  resulting  in  small  peak  height  differ¬ 
ence  particularly  of  the  low-temperature  peaks  which  signif¬ 
icantly  influence  the  thermoluminescence  response  curve 
obtained  by  Gaussian  beam  heating  (see  Case  4  of  this  sec¬ 
tion). 

Case  2:  Polyimide  foil  loaded  with  LiF: Mg,  Ti powder 

Thin  dosimeters  are  desirable  for  applications  in  the 
dosimetry  of  low-energy  0  rays.  We  have  experimented  with 
84-/im-thick  foils  of  polyimide  (Kapton,  DuPont)  which 
were  produced  by  thoroughly  mixing  TLD-100  powder 
(grain  size  20-30  /xm)  with  liquid  polyimide  and  spreading 
this  slurry  uniformly  onto  glass  slides.  After  curing  for  sev¬ 


eral  minutes  at  373  K  and  subsequently  for  1  h  at  573  K,  the 
translucent  yellow  film  was  removed  by  immersion  in  water. 
Correcting  for  the  weight  loss  during  the  polymerization 
process,  the  powder-polyimide  weight  ratio  of  the  dry  foil 
was  about  2:1. 

Pieces  of  2.5  X  2.5  cm2  size  were  mounted  in  metal 
frames  and  exposed  to  y  rays  from  a  “Co  source.  The  ab¬ 
sorbed  dose  was  26  Gy.  Following  a  10-min  preanneal  at  373 
K,  the  thermoluminescence  response  curves  were  measured 
for  a  fixed  laser  power  of  4.78  W  and  various  laser  spot  sizes 
ranging  from  0.039  to  0.098  cm  diameters.  A  family  of 
curves  measured  under  these  conditions  is  depicted  in  Fig.  9. 
Several  different  sites  on  the  same  sample  were  heated  for 
each  laser  spot  size  to  assess  sample  uniformity  and  laser 
stability. 

All  experiments  were  performed  to  test  the  vt'idity  of 
the  theory  presented  in  Sec.  II  B  rather  than  to  fabricate 


FIG.  8.  Measured  and  calculated  (solid  line!  thermoluminescence  response 
curves  obtained  from  the  LiF:Mg,  Ti  (TLD-100)  chip  of  Fig.  7  after  x-ray 
exposure.  The  calculations  were  performed  with  Eq.  |7bl  and  the  first-order 
electron  kinetic  parameters  (frequency  factor,  trap  depth)  determined  by 
McKeever20  for  peaks  3-3  of  this  thermoluminescent  phosphor.  Peak  2  was 
removed  by  10-min  preannealing  at  373  K. 
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FIG.  9  la)  Measured  and  calculated  isolid  line)  thermoluminescence  re¬ 
sponse  curves  obtained  after  x-ray  exposure  of  a  84-^im-thick  polyimide  foil 
loaded  with  LiF  Mg.Ti  (TLD-1001  particles  of  20-30  /im  gram  size  at  a 
polymer-powder  weight  ratio  of  1  2.  The  curves  were  measured  with  a  4,7. 
W  CO-  laser  beam  of  full  widths  ifrom  left  to  right)  of  0.39,  0.62, 0.77,  and 
0.985  mm.  The  calculations  were  performed  with  Eq.  |8)  and  B  =  26 
.  '  K  cdrVs  [see  Eq.  |9|]  which  was  determined  by  the  incandescence  emission 
method  (Appendix  B|  at  this  laser  power.  The  thermal  diffusivity  of  the 
polyimide-LiF  composite  was  found  to  be  a  =  0.0175  cm2/s  by  the  best  fit 
of  the  calculation  to  a  senes  of  measured  curves  obtained  with  vanous  laser 
spot  sizes  but  otherwise  identical  conditions,  lb)  A  set  of  thermolumines¬ 
cence  response  curves  obtained  with  an  identical  sample  as  the  one  used  m 
(al  except  for  longer  heating  times,  which  were  increased  to  demonstrate  the 
occurrence  of  incandescence  emission.  Each  curve  drops  off  sharply  upon 
turning  off  the  laser  beam.  Laser  power:  4.75  W.  Laser  beam  diameter  2w 
Immi  Ifrom  left  to  nght):  0.899,  0.942, 0.985,  1.024,  1.063,  and  1.102. 


viable  dosimeters.  For  example,  Kapton  absorbs  quite 
strongly  the  TLD-100  thermoluminescent  emission  cen¬ 
tered  around  400  nm  and,  therefore,  a  foil  of  this  type  is 
rather  insensitive  as  a  dosimeter.  Nevertheless,  important 
information  on  the  energy  transfer  from  a  laser  beam  to  a 
sample  of  this  type  and  the  resulting  thermoluminescence 
response  can  be  gathered  and  eventually  applied  to  Kapton 
foils  loaded  with  TL  powders  other  than  LiF. 

The  TLD-100  loaded  Kapton  foil  does  not  strongly  ab¬ 
sorb  10.6-/zm  photons.  Yet,  direct  measurement  of  its  ab¬ 
sorption  coefficient  is  not  trivial  because  of  the  photon  scat¬ 
tering  off  the  LiF  particles  dispersed  in  it.  Since  the  foil  is 
quite  thin  and  LiF  has  a  rather  small  absorption  coefficient 
for  l0.6-/xm  photons,  one  can  safely  assume  that  the  heat 
deposition  is  independent  of  z  throughout  the  entire  thick¬ 


ness  I.  However,  scattering  introduces  an  element  of  uncer¬ 
tainty  whose  influence  on  the  agreement  between  theory  and 
experiment  can  only  be  estimated  by  the  comparison  of  mea¬ 
sured  and  calculated  thermoluminescence  response  curves. 
This  case  of  relatively  weak  absorption  is  characterized  by 
fiL<  1.  As  a  result,  Eq.  (7a)  becomes  independent  of  z : 
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In  attempting  to  calculate  the  family  of  curves  in  Fig. 
9(a),  the  problem  of  determining  the  thermal  and  optical 
properties  of  the  sample  that  enter  Eq.  (8)  presents  itself.  For 
example,  it  is  not  immediately  possible  to  obtain  the  relevant 
averaged  thermal  conductivity  k  for  such  a  heterogeneous 
mixture  from  the  known  values  of  its  components.  Also,  the 
total  absorbed  power  P  ( 1  —  R  )( 1  —  e  ~  *L )  is  difficult  to  de¬ 
termine  without  a  suitable  calorimeter,  while  the  averaged 
specific  density  p  and  heat  capacity  c  can  be  calculated  in 
principle  from  the  values  of  the  components. 

However,  the  quantity  B  defined  by  the  entire  param¬ 
eter  combination  in  Eq.  (9)  can  be  measured  by  exploiting  the 
normally  undesired  onset  of  incandescent  emission  that  is 
observed  when  the  sample  is  exposed  to  the  laser  beam  sig¬ 
nificantly  beyond  the  occurrence  of  the  glow  curve  maxi¬ 
mum  [see  Fig.  9(b)],  We  have  developed  a  method  for  this 
task  that  is  based  on  Planck’s  law  of  blackbody  emission.  It 
makes  use  of  he  exponential  time  dependence  of  this  emis¬ 
sion  when  it  is  measured  for  a  narrow  wavelength  interval 
whose  center  wavelength,  A  has  to  be  selected  in  the  red  or 
better  yet  in  the  infrared  so  as  to  not  destroy  the  sample  b  ’ 
overheating  it.  Also,  the  experiment  has  to  be  performed 
with  significantly  higher  laser  power  densities  as  compared 
to  the  ones  used  in  the  thermoluminescence  experiments  in 
order  to  assure  that  the  condition  4at<w2  is  fulfilled  so  that 
measurable  incandescent  emission  is  generated  before  the 
onset  of  thermal  diffusion.  Under  this  condition  the  tem¬ 
perature  rise  of  the  heated  spot  is  linear  in  time  as  can  be  seen 
from  Eq.  (8). 

AT  =  T  —  T0  =  Bt  /w2,  (10) 

and  it  can  be  shown  (see  Appendix  B)  that  the  so-measured 
incandescent  emission  intensity  /  obeys  the  relation 

©/'  =  [a/[Bt /w2  +  T0)2](B /w1),  (11) 

where  a  =  hc/Ak  and  T0  is  the  initial  temperature.  The 
quantity  B  can  be  calculated  by  carefully  measuring  the  laser 
spot  size  2 w  and  dl  /I  as  a  function  of  time.  Details  of  this 
method  are  presented  in  Appendix  B. 

For  the  above-mentioned  TLD-  100-loaded  Kapton  foils 
B/P  was  found  to  be  5.53  (K  cm2/J).  With  this  knowledge 
and  the  known  laser  power  the  thermoluminescence  re¬ 
sponse  curves  of  Fig.  9(a)  can  be  calculated  from  Eq.  (8)  and 
the  computational  procedure  for  glow  curves  obtained  with 
locally  non-uniform  heating  as  described  in  Sec.  II.  This  was 
done  by  normalizing  a  single  member  of  a  set  of  curves  mea- 
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sured  for  different  laser  beam  diameters  2 w  and  selecting  a 
value  for  the  unknown  thermal  diffusivity  a  so  as  to  obtain 
the  best  fit  for  all  members  of  the  set.  Usually  a  rather  good 
fit  is  possible  for  the  curves  measured  with  large  power  den¬ 
sities  of  the  laser  beam  (e.g.,  in  the  present  case  smaller  spot 
sizes  for  the  fixed  laser  power  of  4.70  W).  Larger  spot  sizes 
tend  to  yield  computed  curves  that  rise  later  than  the  experi¬ 
mental  ones  and  are  steeper  prior  to  reaching  the  maximum. 
Perhaps  the  scattering  of  the  laser  light  from  the  LiF  grains, 
disregarded  in  our  computations,  and  the  associated  distor¬ 
tion  of  the  laser  beam  profile  may  be  a  reason  for  this  discrep¬ 
ancy.  The  value  obtained  for  the  thermal  diffusivity 
a  =  0.0175  cmVs  appears  reasonable  in  light  of  the  values 
for  LiF  (0.02743  cnr/s)  and  Kapton  (0.001  cm’/sl. 

In  closing  this  discussion  of  the  powder-loaded  Kapton 
foils  it  may  be  of  interest  to  mention  that  the  measurement  of 
laser-heated  TL  glow  curves  and  the  developed  computa¬ 
tional  procedure  presented  in  this  work  point  to  a  new  meth¬ 
od  for  the  determination  of  the  thermal  properties  of  other¬ 
wise  difficult  to  measure  thin-foil  TLD  dosimeters. 

Having  completely  determined  the  thermal  and  optical 
properties  of  a  given  dosimeter  of  this  type,  it  is  relatively 
straightforward  to  simulate  its  thermoluminescence  re¬ 
sponse  as  a  function  of  other  experimental  parameters  such 
as  laser  power,  beam  diameter,  foil  thickness,  preannealing 
temperature,  and  duration  and  the  like.  As  an  example,  Fig. 
10  shows  the  dependence  on  the  laser  power  of  the  sample 
from  Fig.  9(a)  for  a  fixed  beam  diameter  2 w  =  0.062  cm. 
Thus,  the  techniques  are  now  available  for  the  optimization 
of  dosimeter  designs  for  the  use  in  laser-heated  thermolu¬ 
minescence  dosimetry  readers.3 

Case  3:  Thin  LiF: Mg,  Ti  phosphor  layer  on  a  glass  slab:  (a) 
Heating  beam  on  layer 

Another  dosimeter  configuration  of  potential  impor¬ 
tance  for  laser-heated  thermoluminescence  dosimetry  is  a 
thin  TL  phosphor  layer  on  a  substrate  that  is  transparent  for 


the  emitted  thermoluminescence  photons.  In  this  arrange¬ 
ment  the  heating  beam  may  be  directed  so  as  to  expose  a 
small  spot  of  the  sample  on  the  side  that  is  covered  with  the 
active  layer  and  the  thermoluminescent  emission  is  moni¬ 
tored  with  a  photon  detector  on  the  opposite  side.  The  sub¬ 
strate  may  be  chosen  according  to  its  absorption  coefficient 
for  the  photons  of  the  heating  beam.  For  example  BaF,  LiF, 
and  sapphire  or  glass  are  materials  having  absorption  coeffi¬ 
cients  increasing  in  that  order.  Assuming  the  phosphor  layer 
is  sufficiently  thin  and  composed  mainly  of  TLD- 100  phos¬ 
phor  powder,  little  energy  is  transferred  from  the  10.6-^m 
beam  using  a  BaF  substrate.  However,  useful  thermolu- 
minescence  response  curves  are  obtained  with  a  single-crys¬ 
talline  LiF  substrate  (see  Fig.  1 1).  The  most  efficient  utiliza¬ 
tion  of  the  laser  beam  energy  is  achieved  by  a  highly 
absorbent  substrate  such  as  common  borosilicate  glass.  The 
thermoluminescence  response  of  such  a  dosimeter  configu¬ 
ration  is  the  rather  complex  result  of  both  direct  and  indirect 
(via  absorption  by  the  glass  and  subsequent  diffusion)  heat¬ 
ing  of  the  phosphor.  Both,  its  optical  properties  at  10.6^mas 
well  as  its  thermal  characteristics,  are  due  to  the  contribu¬ 
tions  of  the  individual  properties  of  both  components  in  an 
unknown  fashion  and  are  for  this  reason  difficult  to  measure. 
However,  by  applying  a  layer  that  is  considerably  thinner 
than  the  substrate  and  considering  the  fact  that  this  layer 
does  only  moderately  absorb  the  laser  photons,  heating  is 
mostly  indirect  via  diffusion  from  the  strongly  absorbing 
glass  surface.  Therefore,  one  can  assume  that  the  thermal 
diffusivity  a  of  the  dosimeter  is  very  close  to  that  of  the  glass 
( a  =  0.0075  cm:  /s).  In  fact,  it  should  be  slightly  higher  than 
that  of  lime  gla"  because  of  the  contribution  to  the  “configu¬ 
ration-averaged  '  a  from  the  LiF  content  of  the  phosphor 
layer.  We  have  performed  all  calculations  (see  below)  of  the 
thermoluminescence  responses  of  thin-layer  dosimeter  con¬ 
figurations  which  use  borosilicate  glass  substrates  with  the 
value  of  a  —  0.01  cmVs  (glass  composition:  64%  silicon 
dioxide,  13%  alkali  oxides,  8%  boron  oxide,  and  15%  other 


Time  (seconds) 


FIG.  10.  Example  of  calculated  dosimeter  per¬ 
formance.  The  sample  is  the  self-supported  po- 
lytmide-powder  foil  from  Fig.  9,  all  optical  and 
thermal  properties  being  identical.  The  calcula¬ 
tions  were  carried  out  for  a  fixed  laser  spot  sue 
at  0.62  mm  and  laser  powers  F of  10,  5, 4, 3,  and 
2  W.  Note  that  for  clarity  these  thermolumines- 
cence  response  curves  are  scaled  by  dividing  the 
TL  intensity  by  a  numerical  factor  equal  to  the 
laser  power  in  W. 
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FIG.  11.  Thermoluminescence  response  curves  measured  from  a  uniform 
LiF:Mg,Ti  layer  ithickness:  3.7  mg/cm:.  20-/im  gram  size,  Dow  Coming 
805  silicone  binder)  on  a  0. 15-cm-thick  nonthermolummescent  LiF  single 
crystal  of  25  x  25  mm2  area.  Laser  power:  P  =  4.95  W.  Laser  spot  size  2icun 
mml:  111  0.326.  (2)  0.327.  |3)  0.330,  |4)  0.336.  |5|  0  340  mm.  To  check  the 
reproducibility  of  the  measurements,  several  thermoluminescence  re¬ 
sponses  from  different  sites  on  the  sample  were  measured  for  otherwise 
identical  conditions.  Repeated  measurements  are  shown  here  only  for  the 
largest  three  of  the  five  different  spot  sizes  chosen  for  this  experiment. 


oxides).  The  only  parameter  combination  not  known  for  this 
type  of  dosimeters  is  B  [see  Eq.  (9)],  which  is  obtained  by  the 
best  fit  of  calculated  and  measured  TL  response  curves. 

We  have  performed  experiments  with  thin  uniform 
TLD-100  powder  layers  in  a  silicone  binder  on  0. 159-mm- 
thick  borosilicate  glass  microscope  cover  slides.  An  example 
of  such  a  family  of  glow  curves  measured  with  a  fixed  laser 
power  and  for  various  laser  spot  sizes  is  presented  in  Fig.  12. 
A  25  x  25  mm2  slide  was  coated  with  a  5  l-/zm-thick  layer  of 
LiF  (TLD-100)  powder  (grain  size  between  30  and  40 ^m)  in 
a  silicone  binder  (cured  weight  ratio  2:1)  and  exposed  to  50- 
kV  x  rays.  After  10-min  preanneal  at  373  K  to  remove  peak 
2, 20  well-separated  areas  on  this  sample  were  exposed  to  la¬ 
ser  heating  pulses  of  4.45  W  power  and  spot  sizes  of  full 
widths  ranging  from  1  to  2  mm.  In  order  to  avoid  the  onset  of 
incandescent  emission  from  the  center  of  the  exposed  spot 
the  exposure  times  were  adjusted  according  to  the  power 
density,  ranging  from  100  to  250  ms.  From  the  best  fit  of  the 
computed  curves  to  these  measurements  B  was  found  to  be 
20  K  cm2/s  for  the  laser  power  stated  above. 

With  this  value  of  B  other  series  of  thermoluminescence 
response  curves  for  this  type  of  dosimeter  can  readily  be  gen¬ 
erated.  For  example  the  influence  of  the  laser  power  P  for 
various  spot  sizes  or  the  thickness  L  of  the  sample  can  be 
studied  by  computer  simulation  in  order  to  arrive  at  an  opti¬ 
mal  dosimeter  design  for  a  given  available  CO,  laser.  An 
example  is  presented  in  Fig.  12(b)  which  demonstrates  rather 
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FIG.  12.  lal  Thermoluminescence  response  curves  measured  with  a  4.45- W 
laser  beam  by  heating  the  phosphor-coated  side  of  a  0. 159-mm-thick  glass 
slab  after  x-ray  exposure.  The  51-/jm-thick  coating  consisted  of  TLD-100 
powder  Igratn  size  30-40  /ami  in  a  Dow  Coming  805  silicone  binder  (2: 1 
weight  ratio).  The  solid  lines  are  curves  computed  from  Eq.  (7b)  and  the 
first-order  electron  kinetic  parameters  determined  by  McKeever0  for 
peaks  3-5  Peak  2  was  removed  by  a  10-min  preanneal  at  373  K.  With  the 
known  thermal  diffusivity  a  =  0.01  cmVs  the  parameter  B  was  determined 
by  fitting  a  family  of  curves  obtained  for  the  following  full  widths  Iw  (mml 
(curves  from  right  to  left):  2.03,  1.83.  1.65,  1.45, 1.24.  respectively,  (bl  Com¬ 
puted  thermoluminescence  response  curves  for  various  laser  powers.  The 
sample  is  identical  to  that  of  la).  The  laser  spot  size  2w  is  1.944  mm  and  the 
laser  powers  (curves  from  left  to  right)  are  4.895. 4.45,  and  4.005  W,  respec¬ 
tively. 


large  changes  in  the  peak  location  of  the  thermolumines¬ 
cence  response  curves  when  varying  the  laser  power  by  only 
10%,  indicating  the  necessity  to  use  a  highly  stabilized  laser 
beam  for  applications  of  laser  heating  in  thermolumines¬ 
cence  dosimetry. 

Case  4:  Thin  LiF: Mg,  Ti  phosphor  layer  on  a  glass  slab:  (b) 
Heating  beam  on  glass  side  opposite  the  layer 

It  is  of  course  also  possible  to  heat  the  active  thermolu¬ 
minescent  layer  through  the  substrate.  Choosing  again  a  rel¬ 
atively  thin  microscope  cover  slide  interesting  thermolu¬ 
minescence  response  curves  from  a  TLD- 100  phosphor  layer 
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were  obtained.  Compared  to  those  measured  by  frontal  heat¬ 
ing  with  the  same  laser  power  and  beam  diameter,  they  all 
peak  at  a  later  time  because  diffusion  through  the  strongly 
absorbing  substrate  is  now  the  only  mechansim  of  heat  trans¬ 
fer.  At  the  same  time  the  effected  spot  size  is  considerably 
larger. 

A  series  of  curves  measured  with  a  laser  power  of  4. 18 
W  and  spot  sizes  ranging  from  1.65  to  2. 13  mm  full  width  is 
shown  in  Fig.  1 3  for  the  same  sample  used  in  the  experiments 
of  Fig.  12.  The  theoretical  curves  were  again  computed  from 
Eq.  1 7b)  with  an  experimentally  determined  value  for  the 
parameter  combination  B  =  16.75  cmVs.  The  best  fit  was 
obtained  with  a  thermal  diffusivity  of  a  =  (7.8  ±0.5)x  10" 2 
cm2/s.  The  calculated  dependence  of  the  thermolumines¬ 
cence  response  on  the  laser  power  for  a  spot  size  of  1 .65  mm 
is  plotted  in  Fig.  14  together  with  curves  calculated  for  a 
fixed  power  and  spot  size  and  different  sample  thickness  L 

We  have  used  the  same  sample  in  an  attempt  to  study 
the  effect  of  preannealing  on  the  thermoluminescence  re¬ 
sponse  curve.  These  experiments  are  not  just  a  trivial  repeti¬ 
tion  of  what  has  been  done  before  by  many  authors  to  find 
the  appropriate  treatment  for  removal  of  rapidly  fading  low- 
temperature  glow  peaks  prior  to  the  dose  measurement.  Be¬ 
cause  the  temperature  distribution  is  nonuniform  when  a 
semi-infinite  slab  is  heated  by  a  localized  beam  of  Gaussian 
profile,  uniform  preannealing  effects  the  response  curve 
much  more  severely  than  the  conventional  glow  curve.  This 
becomes  obvious  by  careful  examination  of  the  photographs 
of  Fig.  2  which  may  be  considered  as  a  “coded”  temperature 
distribution  in  the  emitting  layer  for  the  various  times  after 
onset  of  the  laser  exposure.  Imagine  ringlike  surface  areas 
between  two  adjacent  isotherms  of  a  given  small  tempera- 


FIG.  13.  Thermoluminescence  response  curves  of  the  sample  from  Fig.  12, 
however,  the  laser  heating  beam  impinges  onto  the  side  opposite  to  the  phos¬ 
phor-coated  one.  Thermal  diffusion  through  the  thin  glass  slab  is  now  re¬ 
quired  to  produce  the  thermoluminescent  emission,  With  the  known  ther¬ 
mal  diffusivity  of  this  sample  [see  Fig.  12(a)]  the  B  parameter  [Eq.  (9)]  was 
determined  to  be  1 6. 5  K  cm  Vs  at  P*  4. 1 8  W  by  the  best  fit  of  computed 
curves  to  the  entire  senes  of  curves  measured  for  different  laser  beam  widths 
Iw.  Again,  peak  2  of  TLD-100  was  removed  by  annealing.  The  solid  lines 
are  the  curves  computed  with  Eq.  ( 7b)  and  on  the  basis  of  first-order  electron 
kinetics  and  the  following  measured  beam  widths  2w  (mmi:  1  65.  1.994, 
2.13,  respectively. 


FIG.  14.  Thermoluminescence  response  curves  of  the  sample  from  Fig.  1 3 
computed  from  Eq.  i7bl  for  various  laser  powers  /’and  total  sample  thick¬ 
nesses  L  and  fixed  spot  size  2 w  =  1 .65  mm.  Again  it  was  assumed  that  peak 
2  of  the  TLD-100  phosphor  was  removed  by  preannealing.  The  solid  curves 
were  calculated  for  L  =0210  mm  and  the  following  laser  powers  ifrom  left 
tonghti:  6.27.  4.18,  and  3.125  W  The  dashed  curves  demonstrate  the  effect 
of  the  total  thickness  L  by  increasing  the  substrate  thickness  only  (fixed 
phosphor  layer  thickness  of  0  051  mml.  The  calculations  were  performed 
with  a  laser  power  P  =  4. 18  W  and  the  following  total  thicknesses  L  (from 
left  to  right):  0. 105,  0  2 10,  0.420,  and  0.630  mm. 

ture  difference.  At  the  perimeter  of  the  emission  pattern, 
which  is  at  a  low  temperature,  the  emission  from  such  a  ring 
is  due  to  the  low  temperature  peaks  of  the  conventional  glow 
curve.  A  ring  of  the  same  temperature  difference,  but  located 
closer  to  the  center,  emits  at  the  same  time  at  a  higher  tem¬ 
perature.  Assuming  for  the  sake  of  simplicity  identical 
brightness  of  these  two  areas,  their  contribution  to  the  total 
emission  (monitored  with  a  detector  and  without  spatial  re¬ 
solution)  is  proportional  to  the  square  of  their  respective 
radii.  Removing  the  low  temperature  glowpeaks  by  uniform 
preanneal  treatment  affects  the  outer  and,  thus,  larger  emit¬ 
ting  area  and  consequently  the  strongest  contributions  to  the 
emission  at  that  time.  Figure  15  is  an  illustration  of  this 
phenomenon.  Three  thermoluminescence  response  curves 
are  shown  measured  with  a  laser  power  of  3.63  W  and  a  spot 
size  of  1.6-mm  1/e  power  diameter.  The  sample  is  the  same 
as  that  used  to  produce  the  curves  of  Fig.  1 3.  A  373  K  prean¬ 
neal  of  only  5  min,  which  is  known  to  remove  the  low-inten¬ 
sity  peak  2  of  the  conventional  glow  curve20  has  a  rather 
large  effect  on  the  thermoluminescence  response  curve  mea¬ 
sured  with  Gaussian  laser  beam  heating.  Preannealing  at  the 
same  temperature  for  10  min  has  almost  no  measurable  ef¬ 
fect  on  peaks,  3,  4,  and  5  of  the  conventional  glow  curve 
while  the  difference  can  still  be  detected  in  curve  c  of  Fig.  1 5. 
We  refrain  here  from  reporting  preannealing  measurements 
and  computer  simulations  performed  with  all  the  sample 
types  used  in  our  investigations.  Instead  we  restrict  ourselves 
to  pointing  out  that  in  practical  dosimetry  applications  of 
laser  heating  continuous  dosimeter  layers  with  Gaussian 
beam  profiles'9  great  care  has  to  be  taken  to  reproducibly 
preanneal  the  samples.  In  addition,  fading  of  the  information 
on  absorbed  doses  stored  in  the  thermoluminescence  dosi¬ 
meter  will  cause  a  similarly  enhanced  error  as  compared  to 
conventional  contact  heating. 


1636 


J  Appl.  Ptiys  ,  Vol.  56.  No  4.  15  August  1985 


Abtahi  «t  Ml. 


1636 


FIG.  15  Effect  of  preannealing  for  10  min  at  373  K  on  the  thermolumines- 
cence  response  curves  obtained  from  the  sample  of  Fig.  12  by  heating  it 
through  the  glass  with  P  =  3.63  W  and  a  laser  spot  size  2ic  =  1.6  mm  The 
first  curve  ion  the  leftl  is  obtained  without  annealing  immediately  after  ex¬ 
posure  of  the  sample  to  50-keV  x  ray  for  4  min.  The  other  two  curves  repre¬ 
sent  the  TL  response  after  5-  and  10-mm  preannealing  at  373  K,  respective¬ 
ly. 


IV.  CONCLUSIONS 

We  have  presented  the  theory  of  thermoluminescence 
response  curves  generated  with  Gaussian  laser  beam  profiles 
for  a  number  of  practical  and  interesting  dosimeter  configu¬ 
rations,  and  we  have  verified  our  theoretical  approach  by 
direct  comparison  with  experiment.  We  have  also  shown 
that  by  computer  simulation  one  can  assess  the  effects  of 
such  design  parameters  as  layer  thickness,  laser  power  and 
beam  size,  preannealing  temperature  and  duration,  sub¬ 
strate  material  and  its  thickness. 

The  extremely  fast  heating  rates  ( >  104  K/s)  which  can 
be  obtained  by  laser  heating  versus  that  of  conventional  con¬ 
tact  heating  in  solid-state  dosimetry  ( —  1  K/s)  and  the  re¬ 
duced  amount  of  phosphor  required  in  our  thin  film  configu¬ 
ration  indicate  quite  clearly  that  the  techniques  which  we 
have  developed  and  presented  herein  while  of  interest  in 
their  own  right  will  be  of  significant  practical  importance. 


method  is  called  the  Kirchhoff  transformation.32  The  new 
variable  of  transformation  U  is  defined  as 


U 


=  fr*<n 

Jr0  k0 


dT' 


or  by 


^  =  k{T)/k0, 


IA2) 


where  Tis  a  reference  temperature  and  k0  is  the  value  of  it  I T ) 
at  T0.  Upon  substitution,  Eq.  (Al)  becomes 


V2U{rj,t )  +  1  /kog{rj,t )  =  \/a~U (r,z,/ ), 

dt 


(A3) 


where  a  =  k  ( T)/pc . 

For  the  materials  of  interest  in  thermoluminescence  do¬ 
simetry,  the  thermal  diffusivity  a  =  a[T)  varies  very  little 
with  temperature  and  can  be  assumed  constant.  The  initial 
and  boundary  conditions  for  Eq.  (A3)  can  be  derived  from 
the  original  conditions9  using  Eq.  (A2)  and  the  solution  has 
the  same  form  as  Eq.  (7a). 

For  the  case  of  LiF,35  the  thermal  conductivity  may  be 
fitted  by  a  function  of  the  form 

k[T)  =  k0To/T.  (A4) 

In  our  case,  T0  is  the  room  temperature.  Then  the  Kirchhoff 
transformation  gives 

U=T0ln(T/T0),  (A5) 

and  the  temperature  expression  in  terms  of  U  can  be  extract¬ 
ed  from  Eq.  (A5)  as 

T [rj,t )  =  r0  exp[ £/(r,*, f)/r0],  (A6) 

where  U  \rj.t  is  given  as 


U[rj,t  I  = 


1 1  -  R  )aP 
irkL 
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APPENOIX  A 

When  the  thermal  conductivity  varies  with  tempera¬ 
ture,  the  nonlinear  differential  equation  of  heat  conduction 
(compare  with  Eq.  3 1  can  be  written  as 

V  [k\T)VT]  +  g{rj,i)  =  pcdT/dt,  (Al) 

where  the  thermal  conductivity  k  and  the  thermal  diffusivity 
a  are  assumed  to  be  temperature  dependent;  but  the  heat 
source  function  g\rj,t )  does  not  depend  on  temperature. 

Equation  (Al)  can  be  linearized  by  a  transformation 
which  involves  a  change  of  the  dependent  variable.  This 


APPENDIX  B 

The  parameter  combination  defining  2?  in  Eq.  (9)  of  Sec. 
Ill  B  contains  the  thermal  properties,  the  reflectivity,  and 
the  absorption  coefficient  for  the  TL  sample  as  well  as  the 
power  of  the  laser  heating  beam.  In  principle,  these  may  be 
measured  separately.  However,  that  turns  out  to  be  rather 
difficult  for  the  self-supporting  thermoluminescent  foils,  be¬ 
cause  the  material  is  a  composite  of  a  powder  and  a  polymer. 
Therefore,  we  developed  a  method  to  determine  B  experi¬ 
mentally  so  as  to  be  able  to  calcuate  the  temperature  distri¬ 
bution  T[rj,t )  and  the  resulting  glow  curves  of  these  dosi¬ 
meter  configurations,  whose  total  absorption  of  the  laser 
heating  beam  is  not  too  large  (no  dependence  on  z).  It  will  be 
shown  that  it  is  possible  to  extract  B  from  the  incandescent 
background  emission  of  the  heated  spot  on  the  sample. 

As  described  in  Sec.  Ill,  incandescent  emission  from 
the  center  hot  spot  of  a  sample  exposed  to  the  Gaussian  laser 
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beam  is  always  present  when  TL  glow  curves  are  measured 
and  the  laser  beam  is  not  prematurely  turned  off  [see  Fig. 
9(b)]. 

The  energy  density  for  a  given  wavelength  A  of  the  in¬ 
candescent  blackbody  emission  is  given  by  Planck's  law: 

£/,  =  %nhcA  "5[e*P(^)  -  l]  ~ '  (Bl) 

where  T  is  the  absolute  temperature,  and  h  and  k  have  the 
usual  meaning.  We  tried  to  keep  the  temperature  below  600 
K  in  order  not  to  destroy  the  sample  which  visibly  darkened 
by  only  briefly  heating  it  above  that  limit.  Therefore  we  can 
assume  hc/AkT<\,  and  Eq.  (Bl)  reduces  to 


L\  =  8 rrhcA  5  exp 


(- 


he  \ 


AkT. 


(B2) 


The  total  intensity  /  in  the  spectral  range  of  sensitivity  of  the 
photomultiplier  tube  is  given  by 


|B3) 


where  K [A  )  is  the  quantum  efficiency.  It  turned  out  to  be  a 
linearly  decreasing  function  of  the  wavelength  between  600 
and  650  nm,  the  wavelength  interval  utilized  in  our  experi¬ 
ments.  This  yields 


h  a  expl 


<- 


j!~) 

AkT) 


(B4) 


where  A  is  at  the  center  of  this  sensitivity  region  which  is 
determined  by  the  photomultiplier  and  a  short  wavelength 
cutoff  filter  (see  below). 

In  order  to  eliminate  the  unknown  proportionality  con¬ 
stant  in  Eq.  (B4),  one  can  utilize  the  ratio  of  the  first  deriva¬ 
tive  with  respect  to  time  and  this  expression.  This  ratio  is 
readily  determined  experimentally: 


(a/r:)exp(  —  a/T){dT /dt ) 


(B5) 


and 

=  \a/T1)[dT/dt),  <B6) 

where  a  =  hc/Ak. 


The  time  dependence  of  the  temperature  is  given  by  Eq. 
(8)ofSec.  Ill  B.  We  may  consider  the  simplest  case  which  for 
the  temperature  is  proportional  to  the  time  t.  This  is  valid 
only  for  4af  <  w 2  and  yields 


T  —  T0  a  Bt  /w2. 


(B7) 


Then 


[a/(Bt  /w2  +  r0)’]|2?/ur). 


(B8) 


The  quantity  B  is  calculated  by  determining  this  ratio  experi¬ 
mentally  and  careful  measurement  of  the  1/e  full  width  2u> 
of  the  laser  beam. 


EXPERIMENT  RESULTS 

One  of  the  four  different  dosimeter  configurations  was  a 
polyimide  film  loaded  with  LiF  powder.  Both  materials  of 
this  84-/zm-thick  composite  film  absorb  the  10.6-/zm  heating 


FIG.  16.  Determination  of  B  =  ( 1  —  R  l(  1  —  e  ~  “L  )P / irpcL  [see  Eq.  (9)]  for 
powder-loaded  Kapton  foils  by  the  incandescence  method  described  in  Ap¬ 
pendix  B.  Shown  are  the  rise  of  the  incandescence  emission  /  (solid  line)  and 
its  time-denvativeid/  /dt  \/l  (dashed  line)  as  a  function  of  time  after  exposure 
of  the  foil  ( without  previous  exposure  to  ionizing  radiation  to  eliminate  ther¬ 
moluminescent  emissionl  to  a  4.8-W  laser  beam  of  spot  size  2w  =  0. 1 18  cm 
id/  /dt ).  From  the  ratio  id!  /dt  \/I,  determined  between  zero  and  0.03  s.Bwas 
found  to  be  27  +  3  K  cm;/s  for  this  laser  power. 


photons  only  moderately.  One  can,  therefore,  assume  that 
Eq.  (8)  is  applicable. 

The  parameter  combination  B  of  Eq.  (9)  is  straightfor¬ 
wardly  determined  by  recording  the  time  evolution  of  the 
incandescent  emission  at  a  narrow  wavelength  interval 
around  A.  Using  a  step  filter  having  a  low-wavelength  trans¬ 
mission  cutoff  at  600  nm  (interference  filter  type  6135-650, 
Optics  Technology,  Palo  Alto,  California).  Together  with  an 
EMI  9924B  PVIT  yielded  A  =  625  nm.  The  incandescent 
signal  obtained  by  heating  this  foil  with  P  =  4.8  W  and 
2  =  0. 1 1 8  cm  is  shown  in  Fig.  16  together  with  its  derivative 
as  obtained  by  digital  processing  with  the  Data  Precision 
DP6000  waveform  analyzer.  From  the  ratio  \dl  /dt )//,  again 
calculated  directly  with  the  DP6000,  B  was  determined  to  be 
27  cm2  K/s  for  the  laser  power  stated  above.  The  condition 
4af  <  w2  is  fulfilled  with  the  values  a  =  0.0175  cm2/s  (see 
Sec.  Ill  B),  t  =  0.03  sec,  and  w  =  0.059  cm. 

It  should  be  pointed  out  that  the  above  method  is  not 
applicable  for  the  dosimeter  configuration  described  in  Sec. 
Ill,  cases  1,  3,  and  4  because  here  the  z  dependence  of  the 
temperature  rise  cannot  be  neglected  and  Eq.  (8)  is  not  valid. 
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A  general  solution  of  the  heat  diffusion  equation  is  presented  for  the  case  of  a  semi-infinite 
two-layer  system  that  is  heated  with  a  localized  cw  laser  beam  of  Gaussian  or  uniform  circular 
intensity  profile.  As  an  example,  this  theory  is  applied  to  thin  layers  of  a  thermoluminescent 
material  on  glass  substrates.  The  transient  thermoluminescence  emission  response  is  calculated 
and  compared  with  experimental  results,  illustrating  the  validity  of  the  solutions.  Applications 
to  transient  laser  heating  of  oxide-semiconductor  sandwiches  are  discussed. 


I.  INTRODUCTION 

The  interest  in  heating  semi-infinite  slabs  with  localized 
laser  beams  stems  from  the  application  of  laser  heating  in 
thermoluminescence  dosimetry,  which  recently  has  received 
considerable  attention.1'5  The  early  work  in  this  field  was 
restricted  to  the  use  of  cw  CO;  laser  beams  having  Gaussian 
intensity  profiles.  The  theoretical  approach  taken  for  the  cal¬ 
culations  of  the  thermoluminescence  (TL)  response  of  a 
thin  TL  layer  on  a  glass  substrate  assumed  that  the  TL  layer 
thickness  is  negligibly  small  compared  to  that  of  the  sub¬ 
strate.  All  thermal  properties  of  the  sample  were  assumed  to 
be  due  to  the  substrate  alone.  Only  the  reflectivity  and  ab- 
sorptance,  due  to  modification  of  the  glass  substrate  by  the 
thin  TL  layer,  were  taken  into  account.  The  agreement 
between  expenmen tal  TL  responses  and  calculated  ones  was 
satisfactory  for  a  number  of  layer  configurations  consisting 
ofLiF:Mg,  Ti  powder  (Harshaw/Filtrol  TLD-100)  on  bor- 
osilicate  glass.5 

For  practical  applications  of  laser  stimulated  thermolu¬ 
minescence  in  dosimetry  it  became  apparent,  however,  that 
Gaussian  beam  profiles  impose  unnecessary  restrictions.6 
We  have,  therefore,  developed  the  optical  technology  to  pro¬ 
duce  spatially  uniform  beam  intensity  profiles,7  which  was 
modified  for  uniform  circular  beams,  and  we  subsequently 
attempted  to  calculate  the  thermoluminescence  response 
functions  obtained  with  such  a  beam  on  the  basis  of  the  sim¬ 
ple  one-layer  theory  presented  in  Ref.  5.  It  turned  out  that 
this  theory  did  not  yield  the  characteristic  TL  glow  curves 
measured  under  these  conditions.  This  failure  prompted  us 
to  work  out  the  theory  of  localized  laser  beam  heating  of  a 
two-layer  system.  It  is  the  purpose  of  the  present  paper  to 
present  the  general  solutions  of  the  heat  diffusion  equation 
for  this  case  and  apply  the  theory  to  the  same  type  of  thin- 
layer  thermoluminescence  dosimetry  configuration 
(LiF:Mg,TiTL  powder  on  a  borosilicate  glass  substrate) 
used  for  the  work  described  in  Ref.  5.  For  completeness,  we 
have  also  recalculated  the  transient  heating  and  the  resulting 
TL  response  of  this  two-layer  configuration  for  Gaussian 
heating  beams.  It  is  perhaps  worthwhile  to  point  out  that  the 
solutions  described  in  this  paper  are  of  general  validity  and 


may  well  be  useful  in  other  fields  of  endeavor;  i.e.,  laser  an¬ 
nealing  of  semiconductors. 

II.  SOLUTIONS  TO  THE  HEAT  EQUATIONS 

The  system  we  consider  consists  of  two  layers  of  materi¬ 
al  in  perfect  thermal  contact.  As  an  example,  picture  again  a 
thin  layer  of  a  thermoluminescent  phosphor  on  a  glass  sub¬ 
strate  or  an  oxide  layer  on  a  semiconductor.  Each  layer  is 
characterized  by  its  absorption  coefficient  for  the  laser  pho¬ 
tons  /*, ,  thermal  conductivity  k, ,  thermal  diffusivity  a, ,  and 
thickness  /, .  The  index  i  =  1,2  denotes  the  respective  layer 
and  the  total  thickness  is  L  —  /,  +  l2 -  We  also  assume  that 
the  laser  beam  impinges  along  the  z  axis  perpendicularly  to 
the  x,y  plane  of  the  two-layer  configuration.  The  laser  beam 
profile  is  either  circular  and  uniform  (Fig.  1 )  or  Gaussian. 

The  thermal  diffusion  equation  in  cylindrical  coordi¬ 
nates  (ignoring  the  effect  of  the  finite  speed  of  heat  propaga¬ 
tion8)  is  given  by 

r-'d[rk,dT,(rj,t)/dr]  d  [k,dTfr^t)/dz] 

Jr  +  dz 


FIG.  1  Schematic  representation  of  a  uniform  circular  beam  profile  im¬ 
pinging  on  the  surface  at  t  -  0  of  a  two  layer  system  of  total  thickness  L. 
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We  solve  it  for  the  following  boundary  conditions. 

( t )  The  two-layer  slab  is  infinite  in  radial  direction  and 
r,  (rj,t)  —  ro  as  r  —  oo; 

(2)  no  heat  loss  occurs  at  the  boundaries;  i.e..  dTt / 
dz  =  0  atz  =  0(<  =  l )  and  z  =  L(i  —  2); 

(3)  at  the  interface  U  =  /,),  T,  =  T2  and  kxdTx/ 
dz  =  k^dT/dz. 

The  initial  conditions  are  T,  (/%*,;)  =  T0;  i.e.,  room  tem¬ 
perature,  when  the  laser  is  turned  on  at  t  =  0.  The  tempera¬ 
ture  T,  (rjj),  obtained  from  Eq.  ( 1 ),  is  the  actual  tempera¬ 
ture  increase  above  T0.  Let  AT,  =  T,  {rj,t)  —  Tn. 

Correcting  for  reflection  loss  at  the  front  surface  at 
r  =  0,  the  source  functions  become  for  {tP  is  the 

laser  pulse  duration) 

g,(.%z,r)  =  /u , ( 1  -  RX)(P */rrco2) exp (  - fj. ,2 )/( r.ai ) 

(2) 

and 

g2(rj,t)  =n2{\  -  fl,)(l  -  R2)  (P  */rrco2) 

Xexp[  —  >u ,/j -n2(z  -  lx)]f(r,cj)  ,  (3) 

where  f(r,co)  ~  H(cu  —  r)  for  the  uniform  beam  and 
f{r,co)  =  exp(  —  r/oj 2)  for  the  Gaussian  beam.  P*  is  the 
total  laser  power  ( not  the  power  density )  incident  on  the 
front  surface,  H  the  Heaviside  step  function,  R,  the  reflectiv¬ 
ity  of  layer  i,  and  co  the  radius  of  the  uniform  circular  beam  or 
the  1/e  half-width  for  the  Gaussian  beam. 

We  solve  Eq.  ( 1 )  with  the  initial  and  boundary  condi¬ 
tions  given  above  by  the  Green's  function  technique.1*  The 
general  solution,  written  compactly,  is  then 
_ _ _ _ 1 


X  (  dz’  G,,  (rj,t  /r'(z',r)g,(r',r',r) 

Jo 

+  J*  j  dz'  Ga(rj,t  /rj’,T)gz(r'  j!  ,r)^ 
k,  r-  v  dy, 

- - -exp[  -  (yf  +  i7;)(r  -  r)  ] 

a,  Jo  a, 

;V=  V  *.'■ 

,  =  i  2ar,  cos:(£,) 

where  9,  =  £ ,1,/yfa, ' . 


')] 


(4) 


(5) 


x,(z)  = 


cos(£|2/va^)/cos  (£|l|/\ai)  fort 

cos[|2(L  —  z)/^a2\/cos{^2l2/^a2)  for  i 


From  the  boundary  conditions  at  the  interface  we  have 

KiAxT  =  y:/va7  =  P , 


(6) 
=  1 
=  2. 

(7) 

(8) 


and 

=£[ +/9:(a, -a2)  .  (10) 

The  rise  in  temperature  then  becomes 


A T,(rj,t)  =  ( ( 1  —  Rl)P*/irco2N  ] 


f  dr  p  r  dr'  P 
Jo  Jo  Jo 


P  d&  exp[a,02(t  -  r)  ]Jo(0r) 


XJo(0r')f(r',<o)x,(z)<!xp[  -  |  fU  -  r)  ]  {/z[  expf  -MiM 


x[exp(/z,/,)  +  (^l//ulva,)sin(J:1/1/v"a1)  —  cos(^',/1/v'aT)  ]/ 


[(Mi  +  £f/a,)cos(£,/,A/aj]  +  (1  -  R2)fi\ 

Xexp(  -/z,/,)  [cos(£2/2/va2)  +  (£2//z2Va2)sin(£2/2/va^)  -  exp(  -fi2l2)}/ 

[(Mz  +s'2/a2)cos(|'2/;/va^)]}  .  (11) 


This  general  solution  can  be  simplified  for  all  those  two-layer  structures  for  which  one  of  the  layers  is  much  thinner  than 
the  other  one,  which  should  be  valid  for  the  case  of  laser  annealing.  As  we  will  show  below  in  several  examples,  the  spatio- 
temperal  temperature  distribution  of  such  a  configuration  can  be  calculated  by  assuming  a,  =  a2  =  a.  This  assumption 
entails  that  the  general  form  of  the  solution  is  maintained  in  regards  to  the  deposition  of  the  laser  photons  according  to  the 
absorption  coefficients  ju,  and,zz2,  however,  thermal  diffusion  is  governed  mainly  by  the  thicker  one  of  the  two  layers.  Lateral 
heat  diffusion  in  the  thinner  of  the  two  layers  is  negligible.  Instead,  lateral  heating  occurs  via  diffusion  in  the  thicker  layer  with 
heat  transport  across  the  interface  into  the  thinner  layer.  For  this  case,  the  integration  over  P  and  r  in  Eq.  (11)  is  readily 
earned  out,  yielding 

gp  *  r'  /  i  _  n  * 

A 7"  ( rs.i)  = — -  drQ{r,t-r) I  1  -i - — L  V  exp(  -  /z )  [exp(/z,/, )  -  cos(  z/,.,/, )  )sin( J}ml, )  ] / 

u>~  Jo  \  rrtf  ,n  =  , 

[•V-l  1  +  t?i,/Mi:)cos (rjmli)]  +  H  ~  /?;)exp( 

X[cos {rjml2)  +  (Tjm/fJt2)sin(rjmI2)  -exp(  — /z:/2)  ]/ 

Kd  +  Vm/pl  )cos(i)m/:)  ]  jx,  (2)exp[  -ar)l(t-T)]  ,  (12) 
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with  Q(r,t)  =  P(a)/s2at  ,r/^lat  )  for  uniform  circular 
beam  or  exp[ r/(col  +  4 at)  ]/( 1  +  4at  /co2)  for  Gaussian 
beam  profile,  where 

P(cj,r)  =  <j'  dw'  I0(co'  ,r)t\p[  —  (o'2  4-  r)/ 2] 

is  called  the  circular  coverage  function10  or  the  "p"  func¬ 
tion,  1 1 

SiAa,  =sVv«:  =  Vm  • 

B  =  (  1  —  /?,)[!  —  (1  —  /?:)exp(  -/!,/,  —  /I;/;) 

-  R2  exp(  -/u,/,)  }/vNn  , 

and.  from  the  boundary  conditions,  the  transcendental  equa¬ 
tion 

k ,  tan( /, )  +  k2  tan(  t7m  /, )  =  0 ,  (13) 

where  the  roots  77  m .  m  —  1,2...,  are  derived  in  the  Appendix. 

We  illustrate  the  validity  of  these  results  obtained  for 
both  the  uniform  circular  and  the  Gaussian  beam  by  moni¬ 
toring  the  thermoluminescence  response  of  a  thin  TL  layer 
on  a  glass  substrate.  We  choose  as  an  example  a  35-40-/im- 
thick  LiF:Mg,Ti  layer  on  a  1 50-/1  m-thick  borosilicate  micro¬ 


scope  cover  slide.  The  transient  temperature  distnbution  is 
indirectly  monitored  via  the  transient  thermoluminescence 
response  of  the  thin  active  layer  to  localized  laser  heating, 
measured  with  a  photomultiplier.  We  point  out,  however, 
that  in  this  way  the  spatially  integrated  TL  signal  is  mea¬ 
sured  without  spatial  resolution.  Only  the  temporal  evolu¬ 
tion  of  the  thermoluminescence  signal  is  monitored.  The 
procedure  for  calculating  the  thermoluminescence  transient 
is  based  on  the  well-known  electron  kinetics  for  this  thermo¬ 
luminescent  phosphor  (see  below). 

The  absorption  coefficient  of  LiF  for  10.6-/im  CO,  laser 
photons  is  40  cm'1,  while  that  of  glass  is  much  higher 
(~  104  cm-1).  Consequently,  we  have/x,/,<l  for  the  LiF 
layer  and  1  for  the  glass  substrate.  We  have  locally 
heated  this  structure  with  laser  beams  of  Gaussian  and  uni¬ 
form  circular  profiles,  constant  power  P  *,  and  various  beam 
sizes  from  both  sides,  that  is  the  active  TL  layer  can  either  be 
on  top  at  z  =  0  or  at  the  bottom  at  z  —  L,  and  we  have  moni¬ 
tored  the  transient  thermoluminescence  response  in  both 
cases.  The  transient  temperature  distributions  are  readily 
obtained  from  Eq.  ( 12)  as  follows. 

Case  1 .  Thermoluminescent  layer  on  top  of  the  glass 
substrate  (/z  ,/,<!). 


J 


A =  (BP*/co2)  f  dr  Q(r,t-r)  (\  +  Nq  Y  cos(7)mz)exp[  -  arji,  (/ - r)  ]/ 
Jo  V  m  =  1 


[•VmO  *  77^:/^::)cos(7m/,)]) 

Case  2.  Thermoluminescent  layer  at  the  bottom  of  the  glass  substrate  (//,/,<  1). 

A T:(rj,t)  =  (BP*/co2)  f  drQ(r,t-r)(  1  +  ,V0  V  cos[77m  (L  -  z)  ]  e.xp[  —  arj2m  (f  —  r)  ]/ 
Jo  V  nTt  1 

[.VmO  +  )cos(77„,/,)cos(77m/:)  ]  j  . 


(14) 


05) 


In  Eqs.  (12)  and  (15),  we  have  written 

iV0  =  ( l/a)If,  ,  k,I,  and 

Nm  =  ( l/2a)  £  k>l>  { 1  +  Oin  29,  )/(2 9, )  ]/cos2  9,  , 

i  =  \ 

where  9,  —  rjml,.  The  calculations  were  performed  with  the 
following  values  for  the  thermal  conductivities  of  the  two 
layers  in  Case  1  (reverse  in  Case  2): 

( 1 )  TL  layer:  kx  =  0.0275  W/K  cm. 

(2)  Glass  substrate:  k2  =  0.00942  W/K  cm. 

Both  of  these  values  are  temperature  dependent.  The 
temperature  dependence  of  k2  for  glass  is  given  by 
k2  =  0.009 42  (l+crA7\)  with  er  =  0.00222  K-'.'2  This 
can  be  accounted  for  in  our  calculations  by  the  method  of  the 
Kirchhoff  transformation.5  A  “corrected  temperature" 

Tsi  -  T0  +  cr~'[(\  +  2aMz)'n  -  1] 

is  introduced  with  the  same  a.  The  temperature  dependence 
of  k,  for  LiF  does  not  have  to  be  accounted  for  since  it  is  the 
very  thin  layer  and,  as  we  explained  above,  lateral  heat  diffu¬ 
sion  is  dominated  by  the  thicker  glass  layer. 

The  thermal  diffusivity,  a,  assumed  to  be  identical  for 
both  layers  (see  above),  requires  special  attention.  Accord- 
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ing  to  the  glass  manufacturer,12  it  is  slightly  temperature 
dependent  as  well.  It  varies  from  0.0052  to  0.0057  cm:/s 
between  0  and  100  °C,  respectively,  and  it  is  expected  to  be 
somewhat  higher  at  the  highest  temperature  of  400  °C 
reached  in  our  measurements.  We  choose  not  to  correct  for 
this  small  variation  in  a  with  T.  Instead,  we  determine  an 
average  value  for  the  thermal  diffusivity  by  the  best  fit  of  the 
calculate,  and  the  measured  thermoluminescencc  response 
curves  obtained  with  the  uniform  circular  beam.  It  turns  out 
that  the  high  temperature  decay  of  this  curve  is  very  sensitive 
to  slight  variations  in  a.  The  best  fit  to  the  shape  in  all  cases 
reported  on  below  (including  the  ones  where  a  Gaussian 
beam  was  employed )  was  obtained  with  a  —  0.0064  cm2/s. 
This  value  is  consistent  with  expectation  as  the  peak  of  the 
TL  response  curve  occurs  at  about  250  *C,  which  appears  to 
be  a  good  “average”  temperature  for  which  to  select  a  value 
for  the  thermal  diffusivity. 

For  computational  efficiency,  the  numerical  calcula¬ 
tions  of  Eqs.  ( 14)  and  ( 15)  start  at  the  lower  time  limit  of 
t  =  10-”  s  rather  than  zero.  This  can  be  done  since,  for  laser 
powers  that  are  believed  to  be  practical  in  thermolummes- 
cence  dosimetry  applications,  the  temperature  rise  any- 
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where  in  the  sample  during  the  initial  microsecond  is  always 
much  smaller  than  the  temperature  at  which  any  thermolu¬ 
minescence  signal  can  be  detected. 

The  prefactor  B  =  ( 1  —  /?,)( 1  -  R2)/irN0  in  Eqs. 
(14)  and  (15)  is  given  by  the  optical  and  thermal  properties 
of  the  two-layer  sample,  but  it  cannot  be  easily  calculated 
because  of  the  inhomogeneity  of  the  TL  powder  layer. 
Therefore,  we  measure  it  by  the  procedure  described  in  Ref. 
5.  This  procedure  determines  the  absorptance  and  B  is.  ex¬ 
cept  for  a  constant  factor,  the  fraction  of  the  laser  power 
absorbed  by  the  two-layer  system.  Its  measurement  auto¬ 
matically  includes  the  reflection  on  the  front  surface  and  the 
interface.  However,  the  accuracy  of  these  experiments  is 
only  about  ±  25%.  For  calculation  ofTL  responses,  a  value 
within  this  error  is  chosen  to  fit  experimental  data.  The  roots 
obtained  from  the  boundary  condition  are  derived  in  the  Ap¬ 
pendix. 

The  thermoluminescence  glow  curves  we  consider  are 
commonly  called  peaks  3,  4.  and  5.  They  are  sufficiently  well 
described  by  simple  first-order  kinetics  with  attempt-to-es- 
cape  frequencies  of  4x  10' 3  s_ 7.3  X  10 1 5  s  ~ 1 ,  and  4x  1021 
s" and  trap  depths  of  1.23,  1.54,  and  2. 17  eV,  respective¬ 
ly. 13  Their  relative  intensities  are  adjusted  to  yield  the  actual 
glow  curve  measured  with  conventional  uniform  heating. 
Peaks  1  and  2  are  removed  by  a  10-min  pre-anneal  treatment 
at  100  °C.  The  two-layer  configuration  is  irradiated  with 
“Co  y  rays  to  fill  the  electron  traps  and  exposed  to  a  uniform 
circular  CO,  laser  beam  of  fixed  6.5  W  power  both  from  the 
front  (active  TL  layer  exposed  directly)  and  from  the  back 
(active  layer  heated  via  diffusion  through  the  glass  sub¬ 
strate.  which  is  exposed  to  the  laser  beam).  Laser  heating 
generates  a  time-dependent  transient  spatial  temperature 
distribution  which  results  in  the  TL  emission  response 
shown  in  Figs.  2(a)  and  2(b),  measured  without  spatial  re¬ 
solution  by  a  photomultiplier  detector.  As  stated  above,  the 
one-layer  theory  does  not  yield  satisfactory  agreement  with 
measured  thermoluminescence  response  curves  when  the 
heating  beam  has  a  uniform  circular  intensity  profile.  These 
curves  have  a  pronounced  peak:  however,  due  to  lateral  dif¬ 
fusion,  the  thermoluminescence  intensity  does  not  vanish 
completely,  when  heating  is  continued  much  longer  after 
peak  occurrence.  The  simple  two-layer  theory,  Eqs.  (14) 
and  (15),  calculated  both  with  and  without  the  temperature 
dependence  of  the  thermal  conductivity  taken  into  account 
via  the  Kirchhoff  transformation  technique,  is  compared 
with  the  experimental  results  obtained  for  a  fixed  laser  pow¬ 
er  of  6.5  W  and  a  beam  radius  of  a)  =  0.15  cm  [Figs.  2(a) 
and  2(b)].  The  thermal  diffusivity  chosen  was  again 
a  =  0.0064  cnr/s  (an  “average  temperature  value").  The 
prefactor  B  was  selected  to  yield  the  best  fit  of  the  peak  loca¬ 
tion.  It  is  to  be  compared  with  experimental  values,  deter¬ 
mined  with  the  method  of  Ref.  5,  for  both  front  and  back 
heating  of  B  =  10.7  x  ^5%  K  cm:/J  and  B  =  6.3  ±  20% 
K  cm;/J,  respectively  The  agreement  between  the  mea¬ 
sured  and  the  calculated  thermoluminescence  response 
curves  is  very  good,  confirming  the  validity  of  the  theoretical 
approach  presented  above. 

We  have  carried  out  a  detailed  program  exploring  the 
thermoluminescence  response  curves  obtained  under  var- 


Time  (seconds) 


FIG  2.  (a)  TL  response  curve  measured  with  a  6.5-W  CO.  laser  beam 
taj  =  0.15  cm  for  the  uniform  circular  beam  profile)  by  heating  the  phos¬ 
phor  coating  on  tor  >f  a  glass  substrate.  The  dashed  curve  is  obtained  with 
B  -  7  55  K  cm:  J  =  0.0064  cm;/s.  using  Eq.  (14)  of  the  text.  The  full 
curve  is  the  fit  obu  :  id  with  B  =  9  4  K  cm;/J.  the  same  a.  and  the  Kirch¬ 
hoff  transformation  b)  TL  response  when  the  laser  beam  impinges  onto 
the  side  opposite  to  he  phosphor  coating.  The  dashed  curve  is  the  fit  ob¬ 
tained  with  B  =  6  0  K  cm1/}  at  6.5-W  laser  power,  a  =  0.0064  cm;/s.  using 
Eq.  (  15 )  of  the  text  The  solid  line  is  thfe  fit  obtained  with  B  =  10  K  cm  /J. 
the  same  a.  and  the  Kirchhoff  transformation 


ious  experimental  conditions  for  the  purpose  of  characteriz¬ 
ing  this  type  of  laser  heating  in  dosimetry  applications.  The 
results  will  be  presented  in  a  forthcoming  publication.14  15 

Figures  3(a)  and  3(b)  are  families  of  curves  measured 
for  different  1/e  widths  co  of  the  Gaussian  laser  beam.  The 
previous  one-layer  approximation5  yielded  a  satisfactory  fit 
only  by  adjusting  the  prefactor  B.  With  the  two-layer  heat 
equation,  we  obtain  the  curves  drawn  through  the  expen- 
mental  points  in  Figs.  3(a)  and  3(b)  in  one  case,  without 
correcting  for  the  temperature  dependence  of  the  thermal 
conductivity  and,  for  comparison,  when  the  temperature  de¬ 
pendence  is  taken  into  account  by  the  method  of  the  Kirch¬ 
hoff  transformation  as  descnbed  above.  Again,  the  prefactor 
B  had  to  be  adjusted  slightly  for  fit.  We  note  that  the  two- 
layer  theory  represents  a  rather  marked  improvement  of  the 
previous  one-layer  approximation. 

In  conclusion,  we  have  presented  a  general  theory  of 
heating  two-layer  systems  with  localized  laser  beams  of 
Gaussian  and  uniform  circular  intensity  profiles,  and 
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FIG  3.  (a)  TL  response  curves  from  Fig  12.  Ref  5  recalculated  using  a 
4.45-W  CO.  laser  beam  with  a  Gaussian  intensity  protile  b>  heating  the 
phosphor-coated  side  on  a  glass  substrate  The  families  of  curves  are  the  tit 
obtained  with  B  =  5  17  K  cm  J  with  the  Kirchhoff  transformation  (solid 
lines  I  and  B  =  4.05  Kcm;  J  without  the  KirchhofT  transformation 
(  dashed  lines)  using  a  =  0  0064  cm;A.  the  same  as  Fig  2  for  the  following 
full  widths  at  1/e  of  the  peak  intensity  of  the  laser  beam  2w  (mm):  2.03. 
1  83.  1  65.  1  45.  and  1  24  (curves  from  left  to  right,  respectively)  (b)  TL 
response  curves  from  Fig  13  Ref  5  recomputed  with  a  4.18-W  Gaussian 
laser  beam  by  heating  TL  coating  through  the  glass  substrate.  The  dashed 
curves  are  obtained  with  B  =  3.35  K  cnr/J  and  the  solid  lines  are  curves 
obtained  with  B  =  4  3  K  cm  /J.  using  the  KirchhofT  transformation  and 
a  =  0  0064cm;/s  for  the  following  measured  beam  widths  2 ai  l  mm ).  2.13. 
1  964.  and  1  65  (curves  from  left  to  right,  respectively) 

checked  its  validity  by  comparing  calculated  and  measured 
thermoluminescence  response  curves  obtained  with  thin- 
layer  TL  phosphor  layer  on  glass  substrates.  The  theory  is 
well  suited  for  all  systems  that  consist  of  a  layer  of  a  given 
material  in  perfect  contact  with  a  substrate  of  different  opti¬ 
cal  and  thermal  properties.  In  particular,  we  expect  the  the¬ 
ory  described  in  Eqs.  ( 1 )-( 1 1 )  to  be  applicable  to  local  heat¬ 
ing  (laser  annealing)  of  thin  oxide  or  epitaxial  layers  on 
semiconductor  or  insulator  substrates. 
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APPENDIX:  CALCULATION  OF  PARAMETER  i\m 

In  this  work  we  have  a  poorly  absorbing  phosphor  coat¬ 
ing  / 1  =  0.2  L  on  a  highly  absorbing  glass  slide.  l:  =  0.8  L 
where  the  total  thickness  L  =  0.018  75  cm.  In  this  case,  m 
distinct  roots  rjm  can  be  obtained  from  the  transcendental 
equation  tan(0)  -+-  k:  tan(4$)  =  0  with  6  =  0.2rj,„L. 
With  F=5-/L.  '  /j  =  3  » 2/k.  k  =  kl/k~. 

D  -  =  {h  4-  [h1  -  (1  4/Ar )  1  ■  and 

D  ~  =  {h  —  [h  -r-  ( 1  4-  4 /k)  ] 1  :}'  the  m  roots  greater 
than  r)0  =  0  are  then  given  by 

m  “ 1  tan ~'(D  ~  )  +  (m  —  1  )/5]F 

for  m  =  1.6,1 1.  .. 

=  [m-1  tan ~'(D  ~ )  +  (m  —  2)/5]/r 
for  m  =  2.7,12.... 

=  [1  -  jr"'  tan-'(£>  ~)  +  (m  -  3)/5]F 
for  m  =  3,8,13,... 

=  [1  —  it ~ 1  tan (D  ~)  +  (m  -4)/5  f 
for  m  =  4.9,14,... 

=  mF/5 

for  m  =  5,10,15 . 

The  root  =  0  yields  the  first  term,  unity,  in  the  brackets  of 
Eq.  (12). 
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The  theory  of  heating  a  semi-infinite  two-layer  system  with  a  localized  CO:  laser  beam  of 
uniform  circular  intensity  profile  is  applied  to  configurations  consisting  of  thin 
thermoluminescent  LiF:Mg,Ti  phosphor  layers  on  borosilicate  glass  substrates.  We  study 
thermoluminescence  because  it  serves  as  a  convenient  monitor  for  the  temperature  distribution 
and  because  of  its  importance  in  developing  the  laser-heating  technique  for  solid-state 
dosimetry.  Experimental  and  calculated  results  are  compared  in  an  attempt  to  completely 
characterize  the  thermoluminescence  response  curves.  The  two-layer  system  is  heated  in  two 
different  modes:  the  laser  beam  impinges  onto  (a)  the  phosphor  layer,  and  (b)  the  glass 
substrate.  We  investigate  in  detail  changes  in  the  thermoluminescence  response  due  to  different 
laser  powers  and  variations  in  the  optical  and  thermal  properties  of  the  samples. 


I.  INTRODUCTION 

In  a  previous  publication'  computational  methods  were 
presented  for  the  characterization  of  the  thermolumines¬ 
cence  (TL)  response  obtained  from  thin  semi-infinite  ther¬ 
moluminescent  phosphor  layers  on  transparent  substrates 
upon  exposure  to  localized  Gaussian  laser  beams.  These  lu¬ 
minescence  responses  (the  thermoluminescence  intensity 
measured  without  spatial  resolution )  are  not  the  familiar  TL 
glow  curves,2  but  rather  exhibit  a  complex  shape  with  limit¬ 
ed  glow  peak  separation,  which  is  the  result  of  the  nonuni¬ 
form  laser  beam  profile  and  the  resulting  nonuniform  tem¬ 
perature  distribution  in  the  sample.  At  any  given  time  after 
onset  of  laser  exposure,  the  lateral  temperature  distribution 
is  nearly  Gaussian  as  well,  modified  somewhat  by  thermal 
diffusion.  The  spatially  resolved  TL  emission  from  the  thin 
layer  resembles  such  a  Gaussian  distribution  up  to  times 
nearly  2.5  times  longer  than  the  appearance  of  the  peak  of 
the  spatially  unresolved  TL  response  curve. 3  This  curve  typi¬ 
cally  exhibits  a  fast  rise  and  a  shallow  peak,  but  never  decays 
completely  because  of  lateral  heat  diffusion.1  The  TL  signal 
disappears  only  when  the  laser  is  turned  off. 

Nevertheless,  laser  heating  of  thin-layer  thermolu¬ 
minescent  dosimeter  configurations  with  focused  Gaussian 
beams  has  been  applied  successfully  to  TL  imaging4  and  to 
rapid  measurement  of  dose  distributions  generated  by  expo¬ 
sure  to  0  rays  behind  attenuating  filters  of  various  thick¬ 
nesses.5 

For  accurate  and  sensitive  personnel  and  environmental 
thermoluminescence  dosimetry,  heating  with  a  laser  beam  of 
uniform  intensity  profile  is  highly  desirable.  As  we  will  dem¬ 
onstrate  below,  localized  laser  exposure  of  a  semi-infinite 
phosphor  layer  on  a  substrate  produces  glow  curves  which 
are,  except  for  a  slowly  decreasing  tail  on  the  high-tempera- 
ture  side  of  the  peak,  nearly  identical  to  the  ones  obtained 
with  conventional  uniform  heating. 

In  this  paper  we  describe  a  TL  laser  heating  apparatus 
that  features  a  circular  beam  of  uniform  intensity  profile 


generated  with  a  technique  reported  in  Refs.  6  and  7.  It  was 
employed  to  study  the  thermoluminescence  response  curves 
from  thin-layer  LiF:Mg,Ti  dosimeters  on  glass  substrates 
under  a  variety  of  experimental  conditions.  The  motivation 
for  this  investigation  was,  however,  not  limited  to  the  estab¬ 
lishment  of  the  operating  conditions  for  efficient  and  eco¬ 
nomical  laser  heating  in  thermoluminescence  dosimetry  and 
imaging  (the  latter  with  greatly  improved  spatial  resolution 
as  compared  to  the  nonuniform  Gaussian  beam).  We  also 
felt  that  the  thermoluminescence  emission  from  a  thin  layer 
on  a  substrate  is  a  very  convenient  monitor  for  the  time  evo¬ 
lution  of  the  temperature  distribution  established  in  the  two- 
layer  system  during  localized  heating  with  a  laser  beam  of 
uniform  intensity  profile.  The  general  theory  of  the  heat 
transfer  from  such  a  beam  to  a  semi-infinite  two-layer  sys¬ 
tem  and  the  subsequent  temperature  diffusion  was  presented 
previously.3  *  In  the  following  we  will  first  briefly  review  the 
basic  aspects  of  this  theory,  describe  the  experimental  facili¬ 
ty  in  some  detail,  and  finally  present  experimental  and  com¬ 
putation  results. 

II.  COMPUTATIONAL  APPROACH 

In  applying  the  theory  of  Ref.  8  to  a  special  system  tha* 
consists  of  a  thin  thermoluminescent  phosphor  layer  on  a 
strongly  absorbing  substrate,  a  number  of  approximations 
can  be  made  which  greatly  reduce  the  computational  effort 
required  for  simulating  thermoluminescence  response 
curves  measured  under  various  different  experimental  con¬ 
ditions.  The  justification  for  these  assumptions  is  discussed 
in  detail  in  Sec.  IV.  Here  we  only  briefly  review  the  impor¬ 
tant  features  of  the  simplified  theory. 

Assume  a  system  that  consists  of  a  thermoluminescent 
phosphor  layer  on  a  glass  substrate  in  perfect  thermal  con¬ 
tact.  Each  layer  is  characterized  by  its  optical  absorption 
coefficient^, ,  thermal  conductivity  k, ,  and  thickness  /,  The 
index  i  =  1.2  denotes  the  respective  layer  and  the  total  thick- 
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ness  is  L  =  /,  -+-  /,.  In  contrast  to  the  general  theory,  we  as¬ 
sume  for  the  special  case  of  a  thin  phosphor  layer  on  a 
thicker  glass  slide  that  the  thermal  diffusivity  a,  is  tempera¬ 
ture  independent  and  identical  for  both  layers.  The  laser 
beam  has  a  fixed  diameter  and  a  uniform  circular  intensity 
profile  which  is  cylindrically  symmetric  around  the  beam 
axis.  It  enters  the  two-layer  system  perpendicular  to  the  x-y 
plane  either  from  the  side  of  the  phosphor  layer  or  the  glass 
substrate.  The  absorptance  (absorbed  fraction  of  the  inci¬ 


dent  laser  power)  is  of  course  different  in  both  cases,  and  so 
is  the  prefactor  B  in  the  expression  for  the  temperature  rise 
[  Eq.  ( 1 )  below  ] .  The  thermal  diffusion  equation  and  the 
appropriate  boundary  and  initial  conditions  for  this  problem 
are  stated  in  Ref.  8. 

Accounting  for  first-order  reflection  losses  on  the  sur¬ 
face  of  incidence  and  at  the  layer  interface,  the  spatial  and 
temporal  distribution  of  the  temperature  rise  for  the  two 
layers  is  given  bys 


A =  ( BP*/co '-)  j  drQ{r,t-  r)(l  +  [(1  -  Rx)/irB 


x 

X  £  exp( -/r./.Jfexp^,/,) -cos(J7m/,)  +  (»7m/iu1)sin(77m/1)] 

m  =  1 

X  [iV^d  +  rjl'/fi]  )cos(tj„  //, )  ]  4-  (1  -  Rz)e\  p(  -/i,/,) 

X  [cos(»7m/:)  +  (Tjm/n2)sin(rjmlz)  -  exp(  ] 

X  [iVm  ( 1  +  ril/ni  )cos (T)mlz)  ]  “  1  jx,  (c)exp[  -  arj2m(t  -  r)  ]. 


where  P*  is  the  total  laser  power  (not  the  power  density) 
incident  on  the  front  surface.  R,  the  reflectivity  of  layer  i, 
and  a)  the  radius  of  the  uniform  circular  beam.  We  further 
have 

B=  (1  -*0(1  -  (1  -R,)exp(  -Hili) 

—  R:exp(  -nJx)]/wN0, 

N0  =  (1/a)  M,. 


Nm  =  i-  V  k,l,  [  1  +  sin (20,  )/20,  ]cos2(0, ) 

2  ,-1 

with  0,  =  and  Q(r,t)  =  P(co/^2ai ,  r/^lat  ). 

P(oj,r)  =  a)'da}'I0(a)',r)cxp[  —  ( co '2  +  r*)/2] 

is  called  the  circular  coverage  function,  and 

X,(2)  =  COS{7J  m2  )/  COS(  Tf  m  I  j  ) 


X;(2)  =  COS[r?„,  (Z.  -  2)  ] /COS  (?7 

The  parameters  rjm(m  =  1,2,... )  are  the  roots  of  the  follow¬ 
ing  transcendental  equation  (see  the  Appendix): 

fc,tan(?7m/,)  +/c2tan(J7m/2)  =0,  (2) 

which  arises  from  the  boundary  conditions. 

An  additional  simplification  of  Eq.  ( 1 )  is  possible.  It 
stems  from  the  fact  that  the  material  of  interest  in  our  stud¬ 
ies,  LiF,  is  not  a  strong  absorber  for  the  10.6-/zm  photons  of 
the  laser  heating  beam  (p  =  40 cm"*).  In  contrast  to  this, 
the  absorption  coefficient  of  the  glass  substrate  is  high 
(  ~  104  cm  - 1 )  Consequently,  the  laser  deposits  the  energy 
directly  into  the  interface  between  the  layers  if  it  impinges  on 
the  phosphor  side  of  the  two-layer  system  (case  1  below). 
On  the  other  hand,  it  is  possible  to  direct  the  laser  beam  onto 
the  glass  side  of  the  system  where  it  is  completely  absorbed  at 
the  surface  ( case  2 ) .  Thus  the  active  thermoluminescent  lay¬ 
er  can  be  located  between  2  =  0  and  z  =  /,  (case  1)  or 
between  2  =  /,  and  2  =  L.  The  appropriate  expressions  for 
the  temperature  rise  in  these  cases  are  readily  obtained  from 
Eq.  (1). 


Case  1.  Thermoluminescent  layer  on  top  of  the  glass  substrate  (^t ,/,  <  1 ): 

A Tt(rj,t)  =  ( BP*/co 2)  drQ(r,t  -r)^l  +  JV0  £  cos(?7m2)exp[  -  arji  (t  -  r)  ]/[Nm  ( 1  +  17*. //if  ]  j. 


Case  2.  Thermoluminescent  layer  at  the  bottom  of  the  glass  substrate  1 ); 

A7\(r,2,f)  =  (BP*/eo2)  f  dr  Q(r,t  —  t)(\  +  N0  Y  cos[7?m  (Z.  —  z)  ] 

JO  V  I 

Xexp[  —  —  r)  ]/[Nm  ( 1  +  )cos(rfmli)cos(r)ml2)  ]  .  (4) 

The  actual  computer  calculations  of  the  thermoluminescence  response  are  based  on  the  well-known  first-order  electron 
kinetics  of  the  LiF:Mg,Ti  phosphor  (TLD- 100,  Harshaw/Filtrol )  and  require  knowledge  of  the  heating  rates  in  each  volume 
element  Inrdrdz  given  as  follows: 


739 


J.  Appl.  Phys.,  Vol  61,  No.  2. 15  January  1907 


Kelly,  ADtahi,  and  Braunhch 


739 


28 


Case  l  (/*,/,<  1): 

dTy(rj,t)/dt  =  (BP*/co2)Q(.r,t)(\  +  N0  £  cos(77„,z)exp(  -  ari1mt)/[Nm  ( 1  +  rj2mn\  )cos(r}m!l)  ]^.  (5) 

Case  2  (fi2l2  <  1 ) : 

dT2(rj,t)/dt=  (BP*/a>2)Q(r,t)^\  +  N0  cos[vm  (L  -  z)  ]exp(  -  ar)2mt)/[Nm(\  +  tj2„/m]  )cos(^/,  )cos(7m/2)  ]| 

(6) 


The  luminescence  glow  peaks  of  the  LiF:Mg,Ti  considered 
are  commonly  called  peaks  3,  4,  and  5.  They  are  well  de¬ 
scribed  by  first-order  kinetics  with  attempt-to-scape  fre¬ 
quencies  of  4x  10‘\  7.3  X  10‘5,  and  4x  1021  s_1  and  trap 
depths  of  1.23,  1.54,  and  2.17  eV,  respectively. v 

The  sum  of  the  relative  peak  heights'*  of  the  LiF:Mg,Ti 
glow  curve  is  used  as  a  measure  for  the  total  trapped  electron 
population,  and  the  relative  population  of  each  of  the  three 
electron  trap  levels  is  taken  to  be  proportional  to  its  corre¬ 
sponding  peak  height  with  the  correction  of  luminescence 
efficiency  at  different  heating  rates. 10  TL  peaks  that  occur  at 
temperatures  lower  than  those  of  peaks  3-5  are  eliminated 
by  a  preannealing  for  10  min  at  100  °C. 

The  three  first-order  electron-kinetic  differential  equa¬ 
tions  for  peaks  3-5  ( Ref.  9 )  are  solved  with  these  heating 
rates  by  the  Runge-Kutta  method  for  each  volume  element 
2irrdrl{  (for  case  1)  or  2 irrdrL  (for  case  2).  Because  the 
phosphor  layer  is  very  thin,  integration  over  the  phosphor 
layer  thickness  is  omitted  (see  Fig.  10 below).  InSec.  IV  this 


assumption  will  be  justified  by  comparison  of  measured  and 
calculated  thermoluminescence  response  curves.  The  nu¬ 
merical  integration  is  sufficiently  precise  with  step  sizes 
Ar  =  1.25<u/50and  A/ =  10“ 3  s  except  for  the  first  millisec¬ 
ond  during  which  A /  =  10  _6  s. 

III.  EXPERIMENTAL  PROCEDURE 

The  experimental  arrangement  is  schematically  shown 
in  Fig.  1.  A  temperature-stabilized  radio-frequency-excited 
cw  CO;  laser  ( Directed  Energy  Inc.  Model  L35-PS)  emits  a 
beam  whose  intensity  profile  has  an  annular  ( “halo" )  cross 
section  of  1  cm  diameter.  The  10.6-^m  laser  beam  is  linearly 
polarized.  The  maximum  output  power  is  40  W.  For  closed- 
loop  control  of  the  laser  beam,  the  rf  power  supply  is  modi¬ 
fied  for  pulse-width  modulation  at  34  kHz  with  a  variable 
duty  cycle  that  determines  the  laser  output  power. 

The  beam  is  guided  to  a  focusing  module  by  two  parallel 
gold-coated  copper  mirrors.  It  enters  the  light-tight  optical 


light  tight  cnamber 


FIG.  1.  Schematic  experimental  arrangement  of  the  laser  thermoluminescence  apparatus. 
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chamber  through  an  antireflection-coated  Ge  window.  The 
focusing  module  deflects  the  beam  down  onto  the  entrance 
of  the  beam  equalizer  channel,6  which  produces  a  uniform 
square  beam  of  3  X  3  mm2  area.  To  simplify  the  theory  and 
the  numerical  calculation,  the  square  profile  of  the  output 
beam  is  converted  to  a  circular  one  with  the  aid  of  a  2.95- 
mm-diam  aperture  at  the  channel  exit.  The  sample  holder  is 
placed  in  close  proximity  to  it  so  that,  at  the  top  of  the  sam¬ 
ple,  the  beam  has  diverged  only  to  about  3  mm  diam. 

A  gold-plated  electromechanical  shutter,  placed  in 
front  of  the  channel,  controls  the  exposure  time  of  the  sam¬ 
ple  to  the  laser  beam.  When  it  is  closed,  the  laser  beam  is 
reflected  into  a  beam  dump. 

To  provide  a  signal  for  closed-loop  control  of  the  laser 
power,  an  antireflection  coated  ZnSe  window  reflects  a  small 
fixed  portion  (3-5  %)  of  the  laser  beam  to  a  piezoelectric 
detector  (Barnes  Engineering  Inc.  Model  350-5  PZT).  A 
beam-chopping  wheel  operated  at  33.5  Hz  by  a  dc  motor  is 
placed  in  front  of  the  detector  in  order  to  prevent  saturation. 
The  output  of  the  detector  is  fed  to  the  electronic  control 
unit,  providing  an  updated  laser  power  measurement  every 
9.3  ms.  In  this  way,  the  laser  output  power  is  kept  constant 
within  2-3  %. 

A  TL  sample  on  a  25  X  25-mm:  substrate  can  be  placed 
in  a  sample  holder  at  a  distance  of  40-50  /im  from  the  end  of 
the  channel.  A  specific  area  of  the  sample,  selected  by  a  pair 
of  coupled  translation  stages,  can  be  heated  by  the  laser 
beam.  Typical  exposure  times  range  from  100  ms  to  10  s.  The 


translation  stages  are  connected  to  two  motor  driven  actu¬ 
ators  which  are  controlled  by  a  motor  drive  system  with  a 
positioning  accuracy  of  up  to  0. 1  fx m. 

As  in  Ref.  1,  we  have  chosen  LiF:Mg,Ti  (TLD- 100  sup¬ 
plied  by  Harshaw/Filtrol)  as  the  thermoluminescent  mate¬ 
rial  to  serve  as  a  monitor  of  the  temperature  evolution  when 
the  two-layer  system  is  exposed  to  the  laser  beam.  All  the 
systems  which  are  studied  to-date  have  consisted  of  a  35-40- 
/im-thick  fine-grain  phosphor  layer  in  a  high-temperature 
binder  on  150-yum-thick  borosilicate  glass  microscope  cover 
slides.  After  exposure  to  “Co  y  rays  (for  the  purpose  of 
filling  the  traps)  and  a  subsequent  10-min  preanneal  treat¬ 
ment  to  eliminate  low-temperature  glow  peaks,'  the  center 
of  these  samples  was  heated  with  a  stable  and  uniform 
3  +  0.07-mm-diam  C02  laser  beam  with  power  output 
between  2.44  W  and  13.5  W.  The  laser  exposure  was  done  in 
two  different  modes,  the  laser  beam  either  striking  the  front 
(phosphor  layer)  or  the  back  (substrate).  Naturally,  the 
prefactor  B  (see  Sec.  II),  is  different  for  these  two  modes  of 
laser  exposure.  With  the  method  described  in  Ref.  1  its  value 
was  determined  to  be  B  =  10.7  +  25%  K  cm2/J  for  the 
phosphor  layer  surface  and  B  —  6.3  ±  20%  K  cm2/J  for  the 
glass  substrate  surface  ( note  that  B  is  defined  slightly  differ¬ 
ently  since  here  it  is  divided  by  the  laser  power  P  * ) . 

The  emitted  light  from  a  TL  sample  is  transmitted  by  a 
light  guide,3  through  a  shutter,  and  a  protective  glass  win¬ 
dow  to  a  photomultiplier  tube  (PMT,  EMI  type  9924B). 
The  shutter  and  glass  plate  prevent  the  PMT  from  accidental 


FIG.  2.  Thermoluminescence  response  curves  measured  and  calculated  with  a  uniform  laser  beam  of  3  mm  diameter  and  9  75  W  power  heating  at  37  5-^m- 
thick  TL  layer  through  a  0  1 5-mm-thick  borosilicate  glass  substrate  ( the  laser  beam  impinges  onto  the  glass  substrate ) .  The  solid  line  is  computed  with  a 
constant  thermal  conductivity  of  the  glass  k  =  0.009  42  W/K  cm  and  the  dashed  line  with  the  method  of  the  Kirchhoff  transformation,  which  takes  into 
account  the  temperature  dependence  of  k  ( Ref.  1 )  With  the  thermal  diffusivity  a  =  0  0064  cm  Vs  (see  Fig.  3).  the  parameter  B  was  determined  to  be  6  4 
K  cm:/J  by  fitting  to  the  peak  of  the  measured  curve.  This  value  agrees  well  with  the  experimentally  determined  one:  B  =  5.7  t  20%  K  cm  :/L  (see  Ref  1 ) 
Prior  to  the  measurement  of  the  TL  response  curve,  the  sample  was  exposed  to  y  rays  from  a  “Co  source  and  preannealed  for  10  min  at  100  'C. 
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FIG  3.  Thermoluminescence  response  curves  measured  by  heating  with  laser  beams  of  fined  3  mm  diam  but  different  powers  impinging  onto  the  glass  side  of 
the  two-layer  system.  The  sample  is  a  37.5yim-thick  LiF:Mg,Ti  on  a  0. 13-mm-thick  borosilicate  glass  substrate.  The  family  of  curves,  calculated  from  Eqs. 
(4)  and  (6),  yields  the  value  a  =  0  0064  cmVs  for  the  thermal  diffustvity  via  the  best  fit  matching  the  shape  after  peak  position.  The  calculations  were 
performed  with  /c,  -  0.009  42  W/K  cm  that  is  characteristic  for  the  system  and  the  following  B  parameters  and  laser  power:  B  —  6.3, 6.4. 6.0.  and  6.3  K  cm'/ 
3  at  P*  =  13.9  75. 6.5,  and  4.875  W,  respectively  Peaks  1  and  2  of  LiF  Mg.Ti  (TLD-100,  Harshaw/Filtrol )  were  removed  by  10  min  preannealing  at  100  *C 
(Ref.  9). 


exposures  to  visible  light  or  the  unimpeded  laser  power,  re-  The  elect  ronic  control  unit  of  the  apparatus  performs 

spectively.  To  maximize  the  collection  efficiency  of  the  emit-  such  functions  as  timing,  controlling,  and  protecting  the  sys- 
ted  thermoluminescence,  one  end  of  it  is  placed  in  close  tern.  When  a  TLD  reader  heating  cycle  is  activated,  the  fol- 
proximity  to  the  heated  area  on  the  sample.  The  current  lowing  sequence  is  performed: 

signal  of  the  PMT  is  measured,  processed,  stored,  and  dis-  The  laser  is  turned  on  and  its  output  power  is  measured 

played  with  the  aid  of  a  transient  digitizer.  and  adjusted  until  the  laser  is  stabilized  at  the  desired  power 


FIG.  4.  Computed  thermoluminescence  re¬ 
sponse  curves  of  the  sample  from  Fig.  3  for  var¬ 
ious  thermal  diffiisivities  a.  This  fit  with  the 
measured  curve  is  obtained  with  a  —  0.0064 
cmVs  and  the  prefactor  J  =  6.0  K  cmVJ  (dot¬ 
ted  curve) .  The  effect  of  different  thermal  diffu- 
sivities  is  demonstrated  by  the  solid  (a  —  0.005 
cmVs)  and  the  dashed  (a  m  0.0078  cmVs) 
lines.  Laser  power  P*  =*  6.5  W. 
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level.  This  process  is  normally  accomplished  in  100-300  ms; 
however,  it  may  require  a  longer  time  if  the  laser  happens  to 
be  in  an  undesirable  operating  condition.  During  this  time, 
the  shutter  remains  closed  and  the  sample  is  not  exposed  to 
the  beam.  After  the  desired  laser  power  is  obtained,  the  shut¬ 
ter  is  open  for  a  set  desired  heating  time.  At  the  end  of  the 
heating  cycle,  the  shutter  closes  and  the  laser  is  turned  off. 

IV.  RESULTS  AND  DISCUSSION 

In  an  initial  set  of  tests  we  established  the  validity  of  a 
number  of  assumptions  and  approximations  made  to  reduce 
the  computational  effort  required  for  the  comparison  of  ex¬ 
periments  and  the  theory.  The  first  of  these  approximations 
concerns  the  method  of  the  Kirchhoff  transformation  de¬ 
scribed  in  Ref.  1.  It  permits  one  to  account  for  the  tempera¬ 
ture  dependence  of  the  thermal  conductivity  of  the  glass 
substrate,  which  we  approximated  as  k2 
=  0.009  42(1  +  0.002  22A7'; )  W/K  cm,  based  on  data  giv¬ 
en  by  the  glass  manufacturer. 1 1  Figure  2  shows  the  TL  re¬ 
sponse  curve  that  was  measured  by  exposing  the  glass  layer 
(case  2:  heating  from  the  back  of  the  sample)  to  a  9.75-W 
laser  beam  together  with  calculations  based  on  Eqs.  (4)  and 
( 6)  and  in  conjunction  with  the  simplifying  assumptions  dis¬ 
cussed  below.  The  Kirchhoff  transformation  yielded  the 
dashed  line,  and  the  solid  line  was  obtained  assuming  con¬ 
stant  k2  =  0.009  42  W/K  cm  and  with  constant  thermal  dif- 
fusivity  a  =  0.0064  cnrs-1.  In  view  of  the  rather  good 
agreement  between  theory  and  experiment  that  is  achieved 
without  accounting  for  the  temperature  dependence  of  k,  we 
decided  to  disregard  the  Kirchhoff  transformation  for  all 
cases  described  in  this  paper.  We  have  checked  this  with  a 
number  of  different  experimental  situations  and  found  it  was 
always  justified  in  applications  of  the  theory  of  Ref.  8  to  thin 
two-layer  thermoluminescence  dosimeter  configurations. 

The  phosphor  layer  is  indirectly  heated  by  diffusion, 
which  is  rather  rapid  in  the  z  direction,  because  its  thermal 
response  time  or  diffusion  time,12 10  =  l 2 /a,  is  estimated  to 
be  about  1  ms.  This  is  much  shorter  than  the  time  required  to 
measurably  change  the  thermoluminescence  intensity  even 
at  the  highest  laser  power  densities  that  appear  practical  for 
this  type  of  heating  in  thermoluminescence  dosimetry.7  As¬ 
suming  the  thermal  diffusivity  of  the  phosphor  layer  is  not 
significantly  greater  than  that  of  the  glass  substrate,  heat 
transport  by  lateral  thermal  diffusion  occurs  mainly  in  the 
thicker  of  the  two  layers  with  essentially  instantaneous  diffu¬ 
sion  in  the  z  direction  through  the  thin  layer.  This  is  of  course 
also  true  for  case  2  ( laser  heating  from  the  back ) .  Therefore, 
we  are  justified  in  assuming  the  thermal  diffusivity  of  both 
layers  to  be  equal,  its  value  being  that  of  glass.3  Actually, 
there  is  a  reason  to  assume  that  the  thermal  diffusivity  of  the 
phosphor  layer  is  rather  small,  perhaps  even  smaller  than 
that  of  glass.  This  layer  is  a  loosely  packed  matrix  of  small 
grams  in  a  poor  heat  conductor  (Dow  Coming  805  sili¬ 
cone  ) . 

The  thermal  diffusivity  a  for  glass  is  sl'ghtly  tempera¬ 
ture  dependent  (varying  from  0.0052  cmVs  at  0’C  to 
0.0057  cmVs  at  100  ‘C,  but  it  is  unknown  for  higher  tem¬ 
peratures),  and  since  we  did  not  choose  to  account  for  this, 
we  had  the  problem  of  selecting  a  proper  “temperature-aver¬ 


aged”  value  that  provided  a  good  approximation  for  all  cal¬ 
culations  presented  below.  We  established  that  value  as 
a  =  0.0064  cmVs  by  fitting  of  families  of  TL  response 
curves  measured  for  front  and  back  exposure  with  several 
different  laser  powers  ( Fig.  3  is  an  example ) .  This  high  value 
appears  justified  when  one  considers  the  fact  that  the  TL 
glow  peaks  we  selected  as  “temperature  indicators"  for  our 
experiments  reach  their  maximum  emission  intensity  at 
about  250  °C.9  Consequently,  the  major  portion  of  the  ther¬ 
moluminescence  emission  and  thus  the  range  of  “tempera¬ 
ture  monitoring”  occurs  well  above  100  "C.  The  high-tem- 
perature  side  of  the  glow  curve  is  very  sensitive  to  the  value 
of  a  as  shown  in  Fig.  4  and  our  choice  is  the  best  match  to 
this  shape. 


Temperature 


-  2u 


Temperature 


FIG.  5  Computer  simulations  for  the  normalized  spatial  temperature  dis¬ 
tributions  r(rj.t)  calculated  for  exposure  of  a  phosphor  layer  on  a  glass 
substrate  to  a  6.5-W  laser  heating  beam  with  a  uniform  circular  intensity 
profile,  ho  =  3  mm.  The  calculations  are  performed  with  B  =  7  55  K  cm1' 
J,  a  =  0.0064  cm:/s,  and  a  total  thickness  L  =  0. 1875  mm  for  various  times 
after  onset  of  laser  exposure:  (a)  t  =  4  ms  and  (b)  t  =  t  *  =  1 15  ms  (the 
peak  intensity  of  the  thermoluminescence  response  occurs  at  t  * ) 
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Yet  another  simplification  of  the  computational  proce¬ 
dure  is  possible  because  the  TL  phosphor  layer  is  so  thin  and 
diffusion  across  it  is  relatively  rapid.  Instead  of  calculating 
the  transient  TL  emission  from  infinitesimally  thin  layers  5z 
and  integrating  over  z  from  z  =  0  to  z  =  /,,  it  suffices  to  cal¬ 
culate  the  TL  emission  either  from  the  surface  layer  at  z  =  0 
or  from  the  interface  layer  at  z  =  /,.  The  total  instantaneous 
emission  is  then  obtained  simply  by  multiplying  with  the 
thickness,  /,,  of  the  active  layer.  In  each  situation  good  agree¬ 
ment  between  theory  and  measurements  is  obtained.3 
Henceforth,  we  arbitrarily  choose  to  take  the  computed  TL 
emission  from  the  interface  at  z  =  /,  to  represent  the  emis¬ 
sion  from  the  entire  layer. 

We  repeat,  these  simplifications  were  introduced  solely 
to  reduce  the  computational  effort.  The  calculations  have 
yielded  results  that  agree  with  measurements  within  the  ex- 
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penmental  errors.  The  more  exact  calculations  pose  no 
problems  in  principle. 

Even  though  we  have  not  yet  measured  in  detail  spatial¬ 
ly  resolved  thermoluminescent  responses  from  a  TL  phos¬ 
phor  layer  on  a  glass  slide  that  is  heated  with  a  localized 
uniform  laser  beam,  we  present  here  some  initial  computa¬ 
tional  results  as  they  help  to  visualize  the  relation  between 
the  spatial  and  temporal  temperature  distribution  and  the 
corresponding  thermoluminescent  response.  Figure  5(a) 
shows  the  temperature  distribution  in  the  phosphor-glass 
interface  at  f  =  4.  ms  after  onset  of  exposure  of  the  phosphor 
layer  to  a  uniform  3-mm-diam  laser  beam  of  6.5  W  power. 
Figure  5(b)  shows  the  temperature  distribution  at 
t  —  t*  =  1 15  ms  ( here  t  *  is  the  time  of  TL  peak  occurrence ) . 
The  distribution  is  considerably  different  due  to  heat  diffu¬ 
sion.  Interestingly,  the  lateral  distribution  of  the  TL  emis- 
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FIG  6.  The  normalized  spatial  distributions  of  the  thermoluminescence  response  calculated  for  the  sample  from  Fig.  6  at  various  times  t  with  respect  to  the 
time,  t  *.  of  the  TL  peak  intensity  occurrence:  (a)  i  =  U *,  (b)  /  =  r *,  (c)  /  =  1  11 *.  and  (d)  i  =  l.Sf  •. 
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FIG.  7.  Computed  time  evolution  of  the  tem¬ 
perature  of  the  sample  and  experimental  condi¬ 
tions  in  Fig.  6  at  various  distances  r  from  the 
laser  beam  axis:  r  =  0,  o>/4,  <u/2,  to,  and  1. 123  aj 
( from  top  to  bottom ) .  ( a  =  0. 1 5  cm  is  the  radi¬ 
us  of  the  laser  beam. ) 
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sion  retains  the  basically  cylindrical  shape  up  to  t  *;  it  simply 
broadens  with  heating  time.  This  is  of  course  a  result  of  the 
initial  rise  of  a  glow  curve  which  is  proportional  to 
exp(  —  E  /kT)  (£  is  the  thermal  activation  energy).  Imme¬ 
diately  after  t  *,  the  TL  intensity  decreases  in  the  central  and 
hottest  part  of  the  exposed  area  and  becomes  essentially  zero 
at  t  —  \.5t  *,  shown  in  Figs.  6(a)— 6(c) .  For  times  r>1.25f  *, 
the  TL  emission  originates  only  from  those  areas  of  the  phos¬ 
phor  layer  that  are  not  directly  exposed  to  the  laser  beam, 
but  are  heated  by  lateral  diffusion. 

For  applications  of  laser  heating  in  thermoluminescence 
dosimetry  it  is  of  interest  to  investigate  the  temperature  rise 


of  the  exposed  phosphor  layer.  Optimum  operating  condi¬ 
tions  are  commonly  believed  to  include  heating  a  thermolu¬ 
minescent  dosimeter  with  a  constant  heating  rate;  i.e.,  a  tem¬ 
perature  rise  linear  in  time.13  Figure  7  shows  computed 
temperature  rise  curves  for  the  same  sample  and  experimen¬ 
tal  conditions  discussed  above.  In  agreement  with  calcula¬ 
tions  by  Sparks,12  the  temperature  in  the  center  of  the  ex¬ 
posed  spot  rises  first  nonlinearly  for  a  few  ms  and  is  linear 
thereafter.  The  heating  rate  is  proportional  to  the  laser  pow¬ 
er.  The  numerical  calculations  indicate  further,  that  the 
heating  rate  is  very  nearly  the  same  for  an  area  centered 
around  the  laser  beam  axis  up  to  a  radius  0.5tu  for  heating 


FIG.  8.  Calculated  thermoluminescence  re¬ 
sponse  curves  of  the  sample  from  Fig.  3  for  a 
small  variation  in  the  radius  of  the  laser  beam. 
The  computations  are  performed  with  P  •  =  6. 5 
W  and  the  following  beam  radii  a  (from  left  to 
right ) :  0. 1485, 0. 1 5,  and  0. 1 5 1 5  cm. 
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FIG.  9.  A  Mt  of  thermoluminescent  re¬ 
sponse  curves  measured  with  s  37  5-^im- 
thick  TL  layer  on  0. 15-mm-thick  boro- 
silicate  glass  substrate  for  fixed  laser 
beam  radius  oi  =  0. 1 5  cm  and  various  la¬ 
ser  powers  P*.  The  laser  beam  impinges 
on  the  phosphor  side  of  the  sample.  To 
check  the  reproducibility  of  the  mea¬ 
surements,  several  thermoluminescence 
responses  from  different  sites  on  the 
sample  were  measured  for  otherwise 
identical  conditions  with  the  following 
laser  powers  (in  W):  1 3.9.75,6  5. 

4.873,3.25,  and  2.44  (from  left  to  right) 


times  well  past  the  peak  of  the  TL  glow  curve  ( compare  with 
Figs.  4-6).  This  is  of  considerable  practical  importance  as  it 
demonstrates  the  usefulness  of  laser  heating  not  only  for  do¬ 
simetry,314  but  also  in  a  number  of  special  areas  which  ex¬ 
ploit  thermoluminescence,  e  g.,  dating  of  archaeological 
artifacts  and  an  authentication. 15  Very  precise  linear  and 
reproducible  heating  is  possible  of  minute  fine-grain  samples 
that  are  attached  with  a  suitable  binder  to  a  substrate  in  the 
center  of  the  uniform  laser  heating  beam.  The  reproducibil¬ 
ity  and  accuracy  of  this  heating  technique  depend,  of  course, 


on  the  stability  and  reproducibility  of  the  laser  beam.  As 
mentioned  in  Sec.  Ill,  it  is  possible  to  control  the  laser  power 
to  within  3%. 

Experimental  errors  arise  from  the  uncertainty  in  laser 
beam  parameters,  the  phosphor  layer  uniformity,  and  the 
reproducibility  with  which  the  two-layer  samples  can  be  pro¬ 
duced.  With  a  stable  laser  power,  only  the  beam  radius  co 
remained  as  a  critical  parameter.  Calculations  indicated  that 
it  very  sensitively  influences  the  TL  glow  curves.  Figure  8 
shows  that  changing  co  by  only  ±  1%  yields  shifts  in  peak 


0.00  0.10  0.20  0.30  0.40 

Time  (seconds) 

FIG.  10.  Measured  and  computed  thermoluminescence  response  curves  of  the  sample  from  Fig.  9,  and  the  same  experimental  conditions.  The  dashed  curves 
are  computed  from  Eqs.  (3)  and  (5)  and  the  first-order  electron  kinetics  parameters  determined  by  McKeever  (Ref.  9)  for  peaks  3-5  of  LiF  Mg.Ti  and  by 
taking  into  account  the  measured  decrease  in  luminous  efficiencies  with  increasing  heating  rates.  Again,  peak  2  was  removed  by  a  10-min  preannealing  at 
100  *C.  With  the  thermal  diffusivity  a  =  0.0064  cm  Vs  and  a  laser  beam  radius  of  o»  »  0. 1 5  cm,  the  parameter  B  was  determined  for  each  curve  by  the  best  fit. 
The  results  are  calculated  with  the  following  luminescence  efficiencies  17  (in  %  of  the  value  obtained  for  the  smallest  beating  rate):  81.8, 81.9, 85.5, 89.6, 95.7, 
and  100  (from  left  to  right).  The  parameters  B  (k  cmVJ)  are:  8.33,  7.93,  7.55,  7.30,  7.17,  and  7,00  (from  left  to  right). 
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occurrences  that  otherwise  are  produced  only  by  variation  of 
the  produced  layer  thicknesses. 

Variations  in  the  thickness,  /,,  of  the  phosphor  layer 
cause  the  TL  response  curves  to  shift  and  change  in  intensity. 
This  is  demonstrated  in  Fig.  9,  which  shows  a  number  of 
curves  that  were  measured  with  different  values  of  Figure 
9  indicates  the  rather  limited  reproducibility  of  the  spread¬ 
ing  method  used  to  fabricate  the  samples.  Despite  consider¬ 
able  effort,  which  involved  the  construction  of  a  special  layer 
spreading  machine  and  careful  control  of  the  phosphor- 
binder  mixture,  this  performance  could  not  be  improved  be¬ 
yond  the  one  reported  here,  and  more  work  on  layer  fabrica¬ 
tion  techniques  is  needed.  However,  the  layer  uniformity 
achieved  was  sufficient  for  a  study  of  the  TL  response  of  a 
number  of  nominally  identical  samples  as  a  function  of  laser 
power.  Typical  measured  and  computed  glow  curves  result¬ 
ing  from  this  study  as  shown  in  Fig.  10  (case  1:  the  laser 
heating  beam  impinges  onto  the  TL  layer).  The  parameter  B 
had  to  be  adjusted  slightly,  but  note  that  the  variation  in  B  is 
within  the  accuracy  of  the  experimental  method  we  used  to 
determine  it.1  It  was  also  necessary  to  account  for  the  fact 
that  the  overall  luminescent  yield  decreases  with  increasing 
heating  rate  for  otherwise  identical  experimental  conditions. 
Details  of  the  dependence  of  the  yield  on  the  heating  rate  will 
be  presented  in  a  forthcoming  publication.310 

V.  CONCLUDING  REMARKS 

We  have  compared  experimental  and  computed  ther¬ 
moluminescence  response  curves  obtained  from  two-layer 
configurations  by  localized  heating  with  a  3-mm-diam  CO: 
laser  beam  of  uniform  circular  intensity  profile  of  power  up 
to  13  W.  The  results  confirm  the  theory  of  Ref.  8  and  show 
that  thin-layer  thermoluminescence  dosimeters  can  be  heat¬ 
ed  with  heating  rates  up  to  5500  °C/s. 

It  is  of  particular  importance  that  the  thermolumines¬ 
cence  glow  curves  measured  under  these  experimental  con¬ 
ditions  have  the  familiar  shape  obtained  with  conventional 
contact  heating. 
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APPENDIX 

In  this  work  we  have  a  poorly  absorbing  phosphor  coat¬ 
ing  /,  =  Q.1L  on  a  highly  absorbing  glass  slide,  =  0.8L 


where  the  total  thickness  L  =  0.018  75  cm.  In  this  case  m 
distinct  roots  can  be  obtained  from  the 

transcendental  equation  k,  tan  0  +  ic2  tan  4#  =  0  with  6 
=  0.2rjmL.  With  F  -  5 ir/L,  h  =  3  +  2/k,  k  =  6,/fcj, 

Z>+  =  {A+[A2— (l+4/*)]1/2},/2, 
and 

d  -  =  {h-[h  +  (\  +  4/k)]'ny,\ 

then  the  m  roots  greater  than  tj0  =  0  are  given  by 
f [7r-'tan-1(Z>  - )  +  (m  -  1  )/5]F 

for  m  =  1,6,1 1,... 
[fl-_ltan~'(Z>  "■ )  +  (m  —  2)/5]F 

for  m  =  2,7,12,... 

[1  -  rr~  'tan- 1  (D  * )  +  (m  -  3)/5]F 

for  m  =  3.8,13,... 

[  1  —  ir~  'tan  -  '(D  ~ )  +  (m  —  4)/5]F 

for  m  =  4,9,14,... 

mF  /  5 

for  m  =  5,10,15,... 

The  root  ?70  =  0  yields  the  first  term,  unity,  in  the  large 
parentheses  of  Eq.  ( 1 ) . 
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The  general  solution  of  the  thermal  diffusion  equation  for  the  case  of  a  semi-infinite  two-layer 
system  that  is  heated  with  a  localized  cw  laser  beam  of  circularly  uniform  intensity  profile  is 
extended  to  times  after  the  pulse  is  switched  off.  An  application  of  that  solution  in  the  field  of 
solid-state  dosimetry  is  illustrated  by  simulating  the  experimental  thermoluminescence 
response  curves  obtained  from  LiF:  ( Mg, Ti) -phosphor-coated  borasilicate  glass,  for  the  special 
case  of  heating  with  a  C02  -laser  beam  that  has  a  uniform  and  square  intensity  profile:  In 
effect,  we  solve  the  diffusion  equation  experimentally. 


The  Green-function  solution  of  the  transient  diffusion 
equation  was  previously  reported' 2  for  the  case  of  a  semi¬ 
infinite  two-layer  system  heated  by  a  localized  laser  beam  of 
Gaussian  or  circularly  uniform  intensity  profile.  The  solu¬ 
tion  was  valid  for  that  time  /  <  tp ,  where  tp  was  the  pulse 
duration.  The  usefulness  of  the  solution  was  illustrated  by  its 
application  to  laser  heating  with  a  cw  laser  beam  in  the  field 
of  thermoluminescence  (TL)  dosimetry,1'5  where  phos¬ 
phor-coated  borasilicate  glass  comprised  the  two-layer  sys¬ 
tem.  The  predicted  and  experimental  response  was  com¬ 
pared  for  a  wide  range  of  experimental  conditions.  Accepted 
literature  values  were  used  for  the  material  parameters  ex¬ 
cept  for  a  parameter  combination  B  (defined  herein),  which 
was  obtained  by  an  independent  measurement.4  In  fitting 
theory  and  experiment,  B,  the  only  adjustable  parameter, 
agreed  with  its  measured  value  to  within  the  experimental 
error.  However,  in  prior  work,  the  TL  response  was  moni¬ 
tored  only  during  the  pulse  duration  tp,  when  the  source 
term  is  “on”  in  the  diffusion  equation.  Herein,  we  consider 
the  response  of  those  same  samples  but  follow  their  response 
for  some  time  after  the  laser  is  turned  off.  The  generation 
term  is  now  "off”  in  the  transient  diffusion  equation.  As 
well,  the  original  solution  was  only  valid  for  laser  heating 
profiles  that  had  either  a  Gaussian  or  a  circularly  uniform 
intensity  profile:  a  circular  spot.  In  this  note,  we  extend  the 
Green-function  solution12  to  times  after  the  laser  is  turned 
off.  We  then  illustrate  our  procedure  for  analyzing  the  TL 
response  obtained  experimentally  using  a  square  heating 
“spot”  in  terms  of  the  heating-rate  solution  for  a  circular 
spot,  yielding,  in  effect,  a  method  to  solve  the  diffusion  equa¬ 
tion  experimentally. 

The  original  solution  to  the  transient  diffusion  equation 
in  cylindrical  coordinates  (rj)  for  a  two-layer  system  of  to¬ 
tal  thickness  L  =  / ,  +  /  2  yielded  a  heating  rate  as  a  function 


*’  Present  address-  The  Singer  Company,  Singer-Link-Miles.  8895  McOaw 
Road.  Columbia,  MD  20145. 


of  time,  t,  which  is  denoted  herein  as  rh  (ra,t).  It  is  valid  for 
the  boundary  conditions  given  as  follows: 

( 1 )  The  two-layer  slab  is  infinite  in  radial  direction  and 
T(rj,t)  —  7*0  asr—oo. 

(2)  No  heat  loss  occurs  at  the  boundaries;  i.e., 
dT, /dz  —  0  at  z  =  0  (/'  =  1)  and  z  —  L  (i  =  2). 

(3)  At  the  interface  (z  =  /,),  T,  =  T2,  and 
k  /dz  =  k1dT1/dz. 

The  initial  conditions  are  T,(rj,t)  =  T0;  i.e.,  room  tem¬ 
perature,  when  the  laser  is  turned  on  at  t  —  0. 

Correcting  for  the  reflection  loss  at  the  front  surface  at 
z  =  0,  the  source  functions  become  for  0  <  /  <  tp 

g,  ( r,z,f )  =yu,  (1  -  R,  )(/>,/rr<u2)exp(  - /i}z)H{to  -  r) 
and 

gj  (rj,t)  =/ul(l  —  /?,)(!  -  R2)(P*/mo2) 

Xexp[  -/x,/,  —  n2(z  -  /,  )]H(co  -  r), 

where  P*  is  the  total  laser  power  (not  the  power  density) 
incident  on  the  front  surface,  H  the  Heaviside  step  function, 
R,  the  reflectivity  and /x,  the  absorption  coefficient  of  layer  i, 
and  o)  the  radius  of  the  uniformly  circular  beam.  The  rate 
rh  ( rj,t)  is  valid  only  for  times  t  <  tp  where  tp  is  the  time  the 
laser  is  "on”.  If  RH  (rj,t)  denotes  the  heating  rate  for  any 
time  t  where  t>tp>  0,  then  it  is  readily  shown  that 

R„(rj,t)  =rh(rj,t)  -  rh  (rj,t  -  tp)tH{t  -  tp ).  (1) 

For  the  case  of  a  thin  phosphor  layer  on  a  thicker  glass 
substrate  heated  through  the  phosphor,  rh  ( rj,t )  was  given 
explicitly  by 

rh  ( rj,t )  =  ( BP*/(ol)Q(r,t ) 

X  ( 1  +  N0  £  cos(tj mz)enp(  -  arj2mt)/ 

[Nm(\ +T}2m/n\)aMiiml\)])  • 

We  further  have 
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1'  =  (1  —  Rx )  [  1  -  (1  -  R2  )exp(  -fj.xlx  -  H'Jz ) 
-Uj  exp(  -//,/,  )]/irN0, 

N0  =  (i/a)  J  k,I,, 

I  -  | 

and 

Nm  =  t  i 1  +  sin(2e< )/2^-  ]cosz(^, ), 

with 

Q(r,t)  =  P(co/\  2at  ,r/v'2af ), 

k,  the  thermal  conductivity  of  the  layers,  and  a,  =  a;  =  a 
their  approximately  equal  thermal  diffusivity. 

P(a>,r)  =  co'  dco'  /0  (tu\r)exp[  —  (o)'2  4-  r2  )/2] 

is  called  the  circular  coverage  function.  The  parameters  i}m 
(m  —  1,2,.  .  . )  are  obtained  from 

tan(0,  )  4-  k2  tan(#2 )  =  0, 

which  arises  from  the  boundary  conditions. 

Under  certain  experimental  conditions,  it  is  possible  to 
obtain  a  fast  and  linear  temperature  rise.  In  fact,  this  ability 
has  been  used  to  measure  the  luminous  efficiency  of  a  num¬ 
ber  of  different  phosphors.’'6  Then  rh  (r\z,r)  is  a  constant 
and  from  Eq.  ( 1 ) ,  R  H  ( ra,t)  =  0  is  obtained  for  /  >  tp  >  0. 
Equation  ( 1 )  was  derived  assuming  that  the  intensity  profile 
of  the  heating  spot  was  circularly  uniform.  This  solution  for 
small  times  is  not  amenable  to  computation.  However,  for 
short  times,  i.e.,  times  less  than  1  /is,  there  is  no  (useful)  TL 
response,  and  a  good  approximation  is  a  zero  heating  rate  for 
this  time  interval.  We  aiso  note  that  a  TL  response  curve 


simulated  with  this  expression  for  the  heating  rate  has  u» 
maximum  at  a  time  dependent  on  B. 

The  TL  response  of  a  dosimeter  conjuration,  with  the 
heating  rate  given  by  Eq.  ( 1 ) ,  is  simulated  with  the  prior 
procedure.2  It  requires  only  that  the  parameter  B,  the  laser 
power,  and  the  cutoff  time  be  chosen  and  that  an  effective 
'»e  used.  The  TL  response  of  a  dosimeter  configuration  heat¬ 
ed  with  a  square  spot  of  side  is  readily  obtained  experi¬ 
mentally.  This  response  lies  between  that  obtained  for  the 
inscribed  circular  spot  of  radius  co2  =  «  and  of  that  for  the 
circumscribed  circular  spot  of  radius  coc  =  v  2 co.  We  ignore, 
for  the  present  note,  the  temperature  dependence  of  the  ther¬ 
mal  conductivities  k, .  It  has  also  been  shown  that  the  z  de¬ 
pendence  is  minimal  for  a  two-layer  system  when,  as  is  the 
case,  one  layer  is  weakly  absorbing  (/i,  ~40cm  ', 
/i,  —  10* cm “  1  ).  This  remains  true  for  times  t>  tp,  so  we 
choose  the  value  z  —  0.5/ , . 

The  experimental  setup  and  sample  preparation 1 :  have 
been  described  fully.  We  have  repeated  previously  reported 
workz  using  a  modified  version  of  the  TL  reading  apparatus: 
the  circular  aperture  was  removed  (see  Fig.  1  in  Ref.  1 )  to 
produce  a  square  “spot”  of  area  3x3  mm2 .  Figure  1  illus¬ 
trates  the  typical  TL  response  for  several  different  powers 
obtained  for  those  same  samples2  but  for  a  variety  of  laser 
cutoff  times.  TL  response  curves  for  t  =  tp  are  very  similar  to 
those  previously  reported.2 

The  TL  response  for  a  power  of  9.75  W  was  obtained 
with  a  cutoff  occurring  near  the  maximum.  A  fit  to  this  data 
is  obtained  with  Eq.  (1)  with  a  circular  radius  of 
tu,  =  0.3 A/rr  cm.  Only  the  parameter  combination  B  was 
adjusted  to  match  the  location  of  the  maximum  response. 
The  normalized  TL  response  obtained  with  a  uniform  inten- 
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FIG.  J.  Thermoluminescence  response 
curves  measured  (data  points)  and  calcu¬ 
lated  (dashed  lines)  with  a  uniform  laser 
beam  of  3  x  3  mm2  profile  with  the  materi¬ 
al  parameters  described  in  Ref.  2.  The  la¬ 
ser  beam  impinges  on  the  phosphor  sur¬ 
face.  The  curves  from  left  to  right  were 
obtained  with  laser  beam  powers  of  9.75, 
6.5.  4.88,  3.25,  and  2.44  W.  and  corre¬ 
sponding  B  values  of  9.78, 9. 1 3. 8.61, 9. 1 3, 
and  8.76  K  cm2  /J.  The  laser  cutoff  times 
are  80,  128,  200,  251,  and  362  ms.  respec¬ 
tively. 
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sity  profile  should  be  independent  of  the  size  of  the  heating 
“area”  for  semi-infinite  media.  Similarly,  one  should  obtain 
the  same  simulated  response  independent  of  the  choice  of  the 
radius  co.  In  this  regard,  we  note,  for  example,  that  the  simu¬ 
lated  response  for  2.44  W,  which  is  shown  in  Fig.  1,  is  also 
obtained  with  co  equal  to  co,  or  to  coc .  However,  new  values 
are  obtained  for  B  (B  ±  2%,  respectively,  with  B  the  value 
for  co  =  co,).  We  can  eliminate  those  values  of  B  outside  the 
limits  of  experimental  error  for  those  values  of  co  greater  than 
coc  or  less  than  co,  and  get  a  measure  of  the  usefulness  of  this 
approximation.  The  measured  value  of  B  was  determined4 
to  be  10.7  ±  25%  K  cm2/J.  The  values  of  B  obtained  by 
fitting  are  within  this  range  (the  best  fit  being  obtained  with 
co  =  cos ). 

As  Fig.  1  illustrates,  it  is  generally  found  that  the  fit  to 
the  actual  response  is  adequate,  although  there  are  times 
when  it  is  poor,  as  shown  for  the  case  of  4.88  W.  It  is  not  clear 
yet  whether  this  is  due  to  a  nonuniform  phosphor  coating  or 
other  factors.1 2  In  general,  the  simulated  response  is  envel¬ 
oped  by  the  experimental  response,  which  indicates  that  the 
temperature  dependence  of  the  thermal  conductivities 
should  be  included  in  the  comparison.2 

We  have  extended  the  previous  solution1,2  of  the  diffu¬ 
sion  equation  to  all  times,  albeit  for  a  specially  considered 
two-layer  system,  and  we  have  illustrated  its  usefulness  by 
applying  it  in  the  field  of  solid-state  dosimetry.  While  we 
have  considered  a  composite  that  was  heated  through  the 


phosphor  layer,  the  above  analysis  is  readily  extended  to 
such  a  composite  heated  through  the  glass  layer.  It  appears 
then  that  laser  heating  of  practical  dosimeter  configurations 
consisting  of  phosphor-coated  glass  are  adequately  de¬ 
scribed  by  Eq.  ( 1 )  above  provided  only  that  the  glass  is  a 
much  better  absorber  than  the  phosphor.  The  idea  of  using 
TL  response  to  solve  the  diffusion  equation  is  novel.  It  is 
limited  in  this  instance  by  experimental  error  and,  because  of 
the  large  uncertainty  in  B ,  the  present  method  is  not  a  very 
sensitive  test.  Perhaps  a  better  composite  system  can  be 
found:  one  in  which  the  temperature  dependence  of  all  the 
material  parameters  is  known.  Finally,  we  note,  that  the  ap¬ 
proximation  used  herein  is  probably  valid  for  any  shape  of 
the  heating  “spot.” 

We  thank  Norm  L.  Henry,  of  the  Physics  Division  of 
NRC,  for  many  helpful  discussions.  This  work  was  spon¬ 
sored,  in  part,  by  the  Office  of  Naval  Research  under  Con¬ 
tract  No.  N00014-82-K-0529. 
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Part  II:  Summary  of  Work  on  CaSO^Sm 


The  final  phase  of  this  research  program  was  concerned  with  the  optical  stimulation 
of  dose-related  luminescence  from  a  traditional  thermoluminescent  dosimetry  phosphor. 
This  was  achieved  with  some  degree  of  success  with  visible  laser  stimulation  of  a  rare-earth 
doped  CaSC>4  phosphor,  CaSO^Sm.  The  radiation-induced  signal  is  also  optically 
erasable. 

The  essential  features  of  the  process  are  that  normally  the  samarium  occurs  in  the 
triply  ionized  (3+)  state,  but  exposure  of  the  phosphor  to  ionizing  radiation  (we  tested  the 
process  with  ultraviolet  and  gamma  radiation)  results  in  a  dose-dependent  partial 
conversion  of  Sm3+  to  Sm2+  (the  samarium  ions  serve  as  electron  traps). 

Rare-earths  exhibit  characteristic  luminescence  (the  basis  of  several  laser  systems) 
which  can  be  selectively  observed  with  spectroscopic  instrumentation.  Furthermore,  the 
luminescence  spectra  of  the  samarium  ions  possess  excitation  bands  such  that  a  tunable 
stimulating  light  source  can  more  or  less  selectively  excite  the  Sm2+  and  Sm3+  species.  We 
employed  a  tunable,  cw  dye  laser  to  stimulate  the  Sm2+  luminescence  spectrum,  and  an 
optical  multichannel  analyzer  to  record  and  analyze  the  emission  for  dosimetric  information. 

One  interesting  feature  of  the  experiment  turned  out  to  be  that  at  the  chosen 
stimulation  wavelength  (573  nm,  chosen  after  studying  the  excitation  spectrum)  the  Sm2+ 
emission  (==  600-775  nm)  is  more  accurately  described  as  prompt  fluorescence!  The  ion  is 
excited  (by  single-photon  absorption)  to  a  high  level,  and  radiative  relaxation  follows 
promptly.  Thus,  the  readout  or  dose  extraction  procedure  does  not  erase  or  significantly 
degrade  the  dose  information:  the  signal  is  more  or  less  permanently  stored  until  an  erase 
procedure  is  carried  out.  The  signal  can  be  erased  (i.e.,  the  Sm2+  can  be  re-ionized  to 
Sm3+)  by  exposure  of  the  phosphor  directly  to  the  output  of  the  cw  argon  laser  (all-lines) 
that  serves  as  the  dye  laser  pump.  A  full  investigation  of  the  temporal  and  power 
dependent  characteristics  was  not  carried  out.  However,  over  95%  of  the  signal  was 


erased  within  5  seconds  of  the  onset  of  exposure  to  the  sample  to  the  erasing  light, 
although  most  of  this  decrease  occurs  within  the  first  0.5  seconds  [see  Fig.  23,  p.  72  of  the 
Appendix  (M.S.  thesis,  K.  W.  Wells)]. 

Surprisingly,  attempts  to  thermally  reset  the  phosphor  resulted  in  the  opposite  effect 
occurring:  rather  than  reducing  the  radiation  induced  signal,  thermal  treatment  actually 
increased  the  Sm2+  emission.  In  a  similar  fashion,  storage  of  the  phosphor  for  several 
weeks  between  irradiation  and  evaluation  also  led  to  an  increase  of  approximately  equal 
magnitude  compared  with  the  thermal  treatment.  The  signal  does,  however,  stabilize  with 
either  heating  or  time.  We  believe  the  occurrence  of  this  phenomenon  is  the  result  of 
electrons  in  shallow  traps  being  transferred  to  Sm3+  ions.  The  fact  that  heating  the 
phosphor  does  not  erase  the  fluorescence  signal  (even  up  to  600°C  for  1  hour)  points  to  the 
difference  between  this  process  and  thermoluminescence.  The  phosphor  exhibited  a  typical 
TL  behavior:  a  single  transient  glow  curve  is  obtained  upon  heating  the  phosphor  through 
300°C.  Heating  the  material  does  not  affect  those  electrons  trapped  at  Sm2+  sites. 

The  optical  erasure  process  is  dependent  on  the  specific  composition  of  the 
phosphor.  Three  cases  were  studied:  CaSO^Sm  without  any  co-dopant,  CaSC>4:Sm  co¬ 
doped  with  500  ppm  sodium  and  a  third  with  1000  ppm  sodium.  For  all  three 
formulations,  the  radiation-induced  signal  was  completely  erasable  with  the  argon  laser. 
However,  we  did  observe  a  regrowth  of  the  Sm2+  fluorescence  in  the  samples  over  a 
matter  of  hours  after  the  erasure.  We  believe  this  to  be  related  to  the  thermally-induced 
signal  increase  discussed  above.  It  does  appear  that  this  regrowth  can  be  controlled.  The 
greatest  signal  regrowth  occurred  in  the  phosphor  that  was  not  co-doped  with  sodium. 
Addition  of  500  ppm  sodium  cut  this  regrowth  by  approximately  a  factor  of  four,  and  the 
1000  ppm  sodium  formulation  reduced  the  regrowth  by  more  than  a  factor  of  ten  below  the 
case  where  no  sodium  was  intentionally  added.  (Refer  to  Fig.  17,  p.  64  of  Appendix) 
Thus,  it  appears  that  approximately  equal  quantities  of  Sm  and  Na  are  required  to  eliminate 
the  regrowth.  We  did  not,  however,  pursue  this  point  further,  and  the  optimum 


formulation  to  eliminate  the  regrowth  is  as  yet  undetermined,  but  we  believe  attainable. 
Apparently,  the  Sm(3+)  dopant,  which  substitutionally  replaces  the  Ca(2+)  in  the  phosphor, 
needs  to  be  compensated  with  a  singly  positive  ion.  This  compensation  will  occur 
naturally,  or  it  can  be  controlled  by  purposely  adding  an  ion  such  as  Na+.  We  believe  that 
in  this  case,  the  Na+  provides  deep,  stable  electron  traps  that  hold  electrons  which  would 
otherwise  be  transferred  to  Sm3+  to  form  Sm2+,  thus  causing  a  signal  regrowth. 

Dosimetrically,  the  experimental  system  was  found  to  be  very  inefficient  for 
collecting  the  radiation-induced  fluorescence  signal,  and  thus  the  apparent  sensitivity  of  this 
process  is  at  a  level  of  about  0.01  Gray  (1  rad).  However,  the  goal  was  not  to  engineer  a 
high  sensitivity  dosimetry  system  but  rather  to  study  a  possible  process  of  optical 
stimulation  for  recovering  dosimetric  information  on  which  to  base  further  work.  The 
system  can  certainly  be  vastly  improved  to  improve  the  dosimetric  performance. 

One  advantage  of  optical  stimulation  (that  is,  as  opposed  to  thermal  stimulation,  as 
in  laser-heated  thermoluminescence)  is  that  there  is  no  diffusion  of  heat.  Thus,  the  region 
stimulable  with  a  visible  laser  can  be  as  small  as  a  few  microns  in  diameter,  although  we 
have  not  yet  shown  that  this  is  practical.  Thus,  in  principal,  the  "pixel"  size  is  small 
enough  to  produce  high  resolution  images,  or  as  an  erasable  optical  data  recording  medium, 
since  the  phosphor  can  be  "written"  on  and  read  out  with  optical  beams.  A  crude  imaging 
experiment  was  performed  and  is  described  in  Chapter  5  of  the  Appendix. 

Thus,  the  final  year  of  this  ONR  sponsored  program  was  successful  in  achieving  its 
aim:  to  show  that  a  commonly  used  thermoluminescence  dosimetry  phosphor  can  be 
stimulated  by  optical  means  to  yield  dosimetric  information  and  also  can  be  returned  to  its 
initial  state  by  an  optical  beam  to  prepare  it  for  reuse.  The  use  of  a  common 
thermoluminescent  material  is  important  for  the  eventual  acceptance  of  optically  stimulable 
luminescence  in  personnel  radiation  dosimetry.  In  addition,  the  utility  of  optical  stimulation 
will  include  high  resolution  radiation  field  and  dose  distribution  mapping  (e.g.,  in 
radiological  phantoms)  and  make  remote-fiber  optic  dosimetry  more  straightforward. 


Further,  since  there  is  no  glow  curve  involved  in  extracting  the  dosimetric  information, 
there  is  no  apparent  need  to  perform  lengthy,  high  temperature  annealing  procedures  to 
prepare  the  phosphor  for  reuse.  For  this  reason,  it  is  feasible  to  encase  the  phosphor  in 
transparent,  hydrogenous,  low  melting  point  plastic  (e.g.,  polyethylene)  to  serve  as  a 
proton  irradiator  for  neutron  dosimetry  (protons  are  knocked  out  of  the  plastic  by  fast 
neutrons  and  directly  excite  the  phosphor).  Since  no  heating  is  necessary,  the  phosphor 
can  be  read  out  while  still  encased.  The  results  of  this  investigation  are  very  promising  for 
future  developments  in  radiation  dosimetry. 
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